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ABSTRACT 



^ ! Magnetically driven astrophysical jets are related to accretion and involve 



toroidal magnetic field pressure inflating poloidal magnetic fleld flux surfaces. 

Examination of particle motion in combined gravitational and magnetic flelds 
'^ ', shows that these astrophysical jet toroidal and poloidal magnetic flelds can be 

powered by the gravitational energy liberated by accreting dust grains that have 

become positively charged by emitting photo-electrons. Because a dust grain 

K^ I experiences magnetic forces after becoming charged, but not before, charging 

f--. I can cause irreversible trapping of the grain so dust accretion is a consequence 

^^ ' of charging. Furthermore, charging causes canonical angular momentum to re- 

1^^ , place mechanical angular momentum as the relevant constant of the motion. The 

O ' resulting effective potential has three distinct classes of accreting particles dis- 

00 ' 

Q I tinguished by canonical angular momentum, namely (i) "cyclotron-orbit", (ii) 

"Speiser-orbit" , and (iii) "zero canonical angular momentum" particles. Elec- 

"k>( I trons and ions are of class (i) but depending on mass and initial orbit inclination, 

H I dust grains can be of any class. Light-weight dust grains develop class (i) orbits 

such that the grains are conflned to nested poloidal flux surfaces, whereas grains 

with a critical weight such that they experience comparable gravitational and 

magnetic forces can develop class (ii) or class (iii) orbits, respectively producing 

poloidal and toroidal fleld dynamos. 
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1. Introduction 

Magnetohydrodynamically driven plasma jets liaving topology an d dynamics analogou s 



to astrophysical jets have b een produced in l abora tory experiments b ylHsu fc Be 



lanlf2002h : 



Bellan et al. I (120051 ) and by iLebedev et al. I (120051 ) ; see discussion by IBlackmanI (120071 ) . The 
feature of azimuthal symmetry, common to both the lab experiments and to real astrophysical 
jets, has important implications for the structure of the magnetic field. This is because an 
azimuthally symmetric magnetic field can be expressed using a cylindrical coordinate system 
{r, 0, z} as 

B=— (V^X V(/)+Un/V(/)) (i: 



— ( W X V0 + /io/V0) 



where the poloidal flux ?/'(r, z, t) is defined by 

ip{r,z,t)= / 2nr'dr'B,{r',z,t) 
Jo 

and the poloidal electric current /(r, z) is defined by 

I{r,z,t)= / 2Tcr'dr'Jz{r',z,t). 



(2) 



(3) 



The definition of /(r, z) is consistent with Ampere's law for the toroidal field, since using 
V0 = 0/r in EqU] gives the toroidal magnetic field to be 



Ba 



2'Kr 



(4) 



Equations[T]l3]describe the magnetic field and electric current distribution of any axisym- 
metric magnetic field. Because astrophysical jets are azimuthally symmetric, their magnetic 
field must be of the form prescribed by EqslT]- [3] and, indeed, it is generally believed that as- 
tro physic a l jets i nvolve large-scale poloid al magnetic fields threading an accretion disk (e.g.. 



see 



Livid (120021 ): iFerreira fc Cassd (120041 )) and in addition, toroidal magnetic fields. Appli- 



cation of Ampere's law to EqiH shows that the poloidal and toroidal currents are respectively 
given by 



'pol 



and 



1 

2^ 



r^V 



V/x V0 



V^) 



'tor 



27r/io 



V0 



(5) 



(6) 



showing that poloidal magnetic fields are produced by a toroidal electric current and toroidal 
magnetic fields are produced by a poloidal current; toroidal vectors are those vectors in the 
(j) direction and poloidal vectors are those in any combination of the r and z directions. 
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Knowledge of the two stream-function quantities /(r, 2;, t) and iplr, z, t) is thus necessary and 
sufficient to determine the complete vector magnetic field and the complete vector current 
density. 

The term 'magnetic axis' has traditionally been assigned different meanings in the re- 
spective contexts of astrophysics and laboratory toroidal magnetic confinement devices (e.g., 
tokamaks, reversed field pinches, or spheromaks). Specifically, a local maximum in r-z space 
of ip{r, z) is called a magnetic axis in the context of toroidal confinement devices whereas the 
z symmetry axis of the magnetic field is called the magnetic axis in the context of astro- 
physics. To avoid confusion, we will call the location of a maximum of ?/'(r, z) the poloidal 
flux magnetic axis. Poloidal magnetic field lines follow level contours of ip and so one can 
envision the projection of the magnetic field in the r-z plane as being like a set of roads, 
each at a different altitude, encircling a mountain peak at a specific r-z location which is the 
poloidal flux magnetic axis (also called an 0-point). Since a toroidal current at infinity is 
not physical, and since the net magnetic fiux enclosed by a circle with infinite radius must 
vanish as field lines cannot go to infinity, ip must vanish at infinity. Furt hermore, mathe 



matical regularity of physical quantities requires ip to vanish on the z axis (jLewis fc Bellan 



I990I ). Thus, a non-trivial ijj can only be finite in the region < r < 00, —00 < z < cx). 
The simplest situation of physical interest is therefore where ip has a single maximum in 
the r-z plane. We will consider this situation, namely a single poloidal field magnetic axis 
with xpjr, z) symmetri c with respect to z. This situation has been previously considered by 



Lovelace et al. I (120021 ) and implies via EqlHlthat a toroidal current circulates in an accretion 



disk to produce the poloidal magnetic field 

Bpoi= — V^xV0. (7) 

ZTT 

An inescapable feature of this topology is that because V^ = at the maximum of tp, i.e., 
at the poloidal field magnetic axis, the poloidal magnetic field has a null on the poloidal field 
magnetic axis. In the z = plane, the poloidal fiux ip thus starts from zero at r = 0, increases 
to a maximum at the poloidal field magnetic axis, and then decays to zero as r — > cxd. 

We define a to be the radius of the poloidal field magnetic axis. In addition, we define 

{Bz) to be the spatially-averaged axial magnetic field linked by the poloidal field magnetic 

axis and ipQ to be the value of the poloidal magnetic fiux at the poloidal field magnetic axis, 

so 

/;rfr2vrr5,(r,0)_ ^0 ,«, 



{B.) 



J^ dr2nr vra^ 



The axial field Bz = (27rr) ^dip/dr reverses sign at r = a and the radial field Br = 
-{27Tr)~^dip/dz reverses sign at z = 0. An analytic representation for a physically realizable 
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generic flux function satisfying all these properties is derived in Appendix \M This generic 
flux function is 

nr,z)= — : :^^ (9) 



8((^ + l)'+(f)^ 



and has the properties that (i) iplr, z) has a maximum value oiipQ aXr = a^ z = 0, (ii) ip r^ r"^ 
for r <^ a and z = 0, (iii) ip ~ r~^ for r '^ a,z and (iv) for r <^ a/2 or r ^ a and for z ^ a 
the contours of ip are identical to the contours of the poloidal flux produced by a current 
loop located at r = a/2, z = 0. This flux function thus encompasses simpler models which 
assume a uniform axial magnetic field Bz] these simpler models would correspond to the 
r,z <^ a region here since in this region ip r^ r"^ which corresponds to having a uniform axial 
magnetic field B^. This flux function could also be used to describe the far-field of a dipole 
by assuming that r,z ^ a. Since any real axial magnetic field must always be generated 
by a toroidal current located at some finite radius, any real situation will have a poloidal 
flux function qualitatively similar to Eq|9l The flux function presc ribed in Eqj9] is similar in 



essence to the flux function used in Fig.l of iLovelace et al. I (120021 ) . 



Figure [T] plots il){r, z) as prescribed by Eql9]and shows that ^(r, z) has its maximum at 
the poloidal field magnetic axis r = a, 2; = 0. This flux function corresponds to a smoothly 
varying toroidal current density prescribed by Eq|6] concentrated in the vicinity of r = a, z = 
0. Since for 2; = and small r, this function has the asymptotic dependence if) ~ 27?/'o ('"/ct) 1 
it corresponds to an approximately uniform axial magnetic field Bz — 27?/'o/vra^ for r,z <^ a. 
The r <^ a inner-region Bz is thus 27 times stronger than the average Bz field between 
and a. The total toroidal current X^ associated with the generic flux function given by Eq|9] 
is calculated in Appendix [B] using the integral form of Ampere's law and found to be 

X, ^ i!*l . (10) 

vra/io 

The laboratory jets involve the mutual interaction between poloidal and toroidal mag- 
netic fields powered by laboratory capacitor banks. The jet acceleration mechanism results 
from the pressure of the toroidal magnetic field inflating flux surfaces associated with the 
poloidal magnetic field. The question arises as to what powers the toroidal and poloidal 
magnetic fields in an actual astrophysical situation. Existing models of astrophysical jets 
are based on the magnetohydrodynamic (MHD) approximation of plasma behavior and typ- 
ically assume (i) the poloidal field is pre-existing and (ii) the toroidal field results from a 
rotating accretion disk twisting up this assumed primordial poloidal field. The purpose of 
this paper is to present an alternate model wherein it is postulated that the toroidal and 
poloidal field result instead from a non-MHD dusty plasma dynamo mechanism that converts 
the gravitational energy of infalling dust grains into an electrical power source that drives 



-5- 



Fig. 1. — Plot of the normalized generic flux function ■?/'(r, z)/iIjo in coordinates normalized to 
the radius of the poloidal field magnetic axis, that is to the radial position of the maximum 
of i/j{r, z) . Contours of iso-surfaces shown on top; these correspond to projection of poloidal 
magnetic field onto r-z plane. 



o 




-6- 



poloidal and toroidal electric currents creating the respective toroidal and poloidal fields. A 
brief o u tline of how infalling charged dust can drive poloidal currents has been presented in 
BellanI ^mk . 



This model obviously requires existence of sufficient infalling dust to provide the jet 
power. Since the dust-to-gas mass ratio in the Interstellar Medium (ISM) is 1%, one might 
be tempted to argue that any jet driven by the proposed dust infall mechanism would be 
limited to having less than 1% of the powe r available fr om infalling gas, a constraint that 
would contradict observations. However, in iBellaru (120081 ) (to be referred to as Paper I), we 
showed that the dust-to-gas mass ratio in a molecular cloud can be substantially enriched 
compared to the ISM value (e.g., the dust to gas mass ratio in a molecular cloud could be 
enriched 20-fold compared to the 1% ISM value). This enrichment occurs because accreting 
dust slows down much more in proportion to its initial velocity than does accreting gas so 
that the density amplification resulting from dust slowing down is much greater than the 
corresponding density amplification of gas. 

The condition for the toroidal magnetic field to infiate the poloidal magnetic field and 
create a jet can be expressed as 

^>A (11) 

where A is a parameter of the order of the inverse characteristic linear dimension in the 
radial direction. The ratio I/ip can be thought of as the ratio of the electric current fiowing 
along a fiux tube to the magnetic fiux content of the fiux tube and is proportional to the 
twist of the rnagnetic field. Equat i on [TTl well-established in spheromak formation physic s 
jBarnes et al. I Il990l : Ijarboel llQQJ ICeddes et al\ llQQsl : iBellanI I2OO0I : IHsu k BellaJboOsh . 
is essentially a statement that jet expansion (i.e., poloidal field infiation) occurs when the 
toroidal magnetic field pressure force ~ B^Ai acting on area Ai exceeds the restraining force 
BIA2 of the poloidal magnetic field 'tension' acting on area A2. Here Ai and A2 are not 
exactly the same because the toroidal and poloidal fields do not act over the same areas. 
The equivalence between Eq. [11] and the condition B"^ > BIA2/A1 is seen by substituting 
/iq/ = 2'KaB^ from Ampere's law and if) ~ B^ira'^ in EqJTTl 

Paper I divided the regions of interest into successively smaller concentric regions and 
considered dust and gas behavior in the outermost regions. Simultaneous gas and dust ac- 
cretion were considered and it was shown that the dust could be considered as a perturbation 
on the gas, so that the gas accretion problem could be solved first without considering dust 
and then the solution of this gas accretion problem could be used as an input for the dust ac- 
cretion problem. Below is a listing showing which regions were considered in Paper I, which 
are considered in this paper, and which will be considered in a future paper; the nominal 
radii scales and star mass are from Table 3 in Paper I: 
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ISM scale (considered in Paper I): The outermost scale is that of the Interstellar Medium 
(ISM). The ISM has a gas density ~10^ m~^, a dust-to-gas mass ratio of 1 percent, 
a gas temperature T^^'^^ ~ 100 K, and is optically thin. The ISM is assumed to be 
spatially uniform and to bound a molecular cloud having radius r^dge ~ 10^ a.u. 

Molecular cloud scale (considered in Paper I): The molecular cloud scale has much higher 
density than the ISM and is characterized by force balance between gas self-gravity and 
gas pressure. The molecular cloud scale is sub-divided into a large, radially non-uniform 
low-density outer region and a small, approximately uniform, high-density inner core 
region. Clouds have a characteristic scale given by the Jeans length rj ~ 1.4 x 10*^ a.u. 
The radial dependence of gas density is provided by the Bonnor-Ebert sphere solution 
which acts as the outer boundary of the Bondi accretion scale. 

Bondi accretion scale (considered in Paper I): The Bondi accretion scale is ~ r^ ~ 4.3 x 10'^ 
a.u. which is sufficiently small that gas self-gravity no longer matters so equilibrium 
is instead obtained by force balance between gas pressure and the gravity of a central 
object assumed to be a star having mass M ~ 0.4Mq The Bondi scale is sub-divided 
into three concentric radial regions: an outermost region where the gas flow is subsonic, 
a critical transition radius at exactly r^ where the flow is sonic, and an innermost region 
where the gas flow is free-falling and supersonic. 

Collisionless dusty plasma scale (considered in this paper): Free- falling dust grains collide 
with each other in one of the above scales and coagulate to form large-radius grains 
which are collisionless and optically thin. The optically thin dust absorbs UV photons 
from the star, photo-emits electrons and becomes electrically charged. The charged 
dust grains are subject to electromagnetic forces in addition to gravity. Motions of 
charged dust grains relative to electrons result in electric currents with as s ociate d 



poloidal and toroidal magnetic fields [see preliminary discussion in iBellanI (120071 )]. 
This region is assumed to have a scale of 10 — 10^ a.u. and corresponds to the scale 
of a, the radius of the poloidal magnetic field axis. It is assumed that a distributed 
toroidal current peaked at a nominal radius a is responsible for producing a poloidal 
field having the generic profile given in Fig{Tl 

Jet scale (to he considered in a future publication) : The electric currents interact with the 
magnetic fields to produce magnetohydrodynamic forces that drive ast rophysical jets 



i n a manner con sistent with E g ITT] an d analogous to that reported in iHsu fc Bellan 



3 



(I2OO2I . I2OO5I ) and iBellan et al. I (120051 ). This region is assumed to have a scale -C 10 
a.u., possibly as small as a few a.u. and will involve a deformation of the generic 
poloidal field profile given in FigH] because of the pressure of toroidal magnetic field 
infiating the poloidal fiux surfaces. 



Outline of model 



We will show how infalling collisionless dust grains can develop special three dimensional 
orbits suitable for sustaining both toroidal and poloidal dynamos. This result is obtained 
by considering Hamiltonian particle dynamics in the combination of the gravitational field 
of a star with mass M and a three-dimensional axisymmetric magnetic fie ld topology c on- 



sistent w it h pre vi ous models of m a gneti cally d riven astroph y sical j ets [e.g.. iLovelacd (Il976l ). 
Li et al. I (I2OOII ). iLovelace et al. I (120021 ). and iLynden-Belll (120031 )]. The reason why dust 
grains develop these special orbits will be shown to be due to charging of dust grains via 
photo-emission of electrons. The analysis involves using Hamiltonian mechanics to gener- 
alize the centrifugal potential so as to include magnetic force, i.e., the Stormer effective 
potential is used. Storm er potentials have been previously used for invest igating auroral 



particles (IStormerlll955l ). electron and ion motion in the magnetosphere (jShebalinl 12004 



Lemaird l2003l ) and most rece ntly, charged dust grain motion in the m agneto-gravitational 



fields of Saturn and Jupiter (JDuUin et al. I l2002l : [Mitchell et al. 1 120031). Stormer poten tials 



are also commonly used to characterize particle orbits in tokamaks (JRome fc Peng||l979l ) and 
St ormer potentials were found to be important in the MHD-driven jet experiment reported 
by JTripathi et al. I (120071 ). We will restrict the analysis to showing how toroidal and poloidal 
field dynamos can be sustained in steady state by these special Hamiltonian particle orbits; 
the much more complicated problem of how a dynamo grows from a seed magnetic field will 
not be addressed here. These special orbits are quite different from conventional cyclotron 
orbits. As reviewed in Appendix [C] a dynamo cannot be sustained by particles executing 
cyclotron orbits because cyclotron orbits and associated drifts are diamagnetic, i.e., create 
magnetic fields that oppose the field in which the particle is orbiting. 

The importance of a Hamiltonian analysis can be appreciated by considering the gedanken 
experiment where the charge to mass ratio of a particle in a combined gravitational-magnetic 
field is assumed to be increased from zero (neutral particle) to that of an electron or ion. The 
particle will thus make a transition from Kepler to cyclot r on orb ital motion. The details of 
how this transition occurs have been examined by JBellani (120071 ) in the context of uniform- 
magnetic-field orbits restricted to a plane. The present paper will address this issue in the 
more general context of three dimensional particle orbits in a spatially non-uniform three 
dimensional magnetic field having dipole-like topology appropriate for an accretion disk; 



simila r dipole topology has been previously invoked for accretion disks by iLovelace et al. 
(I2OO2I ). Our analysis identifies five distinct classes of orbits and shows that the class to 



which a given charged particle belongs depends both on its charge to mass ratio and on 
the circumstances under which the charged particle was created from an initially neutral 
particle. The interaction between the distinct symmetries of the magnetic and gravitational 
fields removes the isotropy of the incident neutral particles existent prior to charging so 
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that the newly formed charged particles separate into groups having qualitatively different 
types of orbits. Some orbits correspond to a simple accretion, some involve accretion and 
production of a dynamo driving toroidal current, and some involve accretion and a dynamo 
driving poloidal current. The type of orbit a charged particle develops depends on both the 
angular momentum and the angle of incidence of the parent neutral particle. 

The paper is organized with the goal of being concise while also realizing that some read- 
ers may not be familiar with the concepts of adiabatic versus non-adiabatic orbits, Speiser 
orbits, Stormer effective potentials, and how conservation of canonical angular momentum 
results in confinement of an adiabatic particle to the vicinity of a poloidal flux surface. 
Rather than reviewing these concepts in an introductory section , they are instead discussed 
in appendices. 



3. Reduction of collisionality due to dust agglomeration 

Paper I showed that dust grains are collisionally decoupled from gas in the ISM and 
then become collisionally coupled to gas in the Bonner-Ebert and Bondi regions of a molec- 
ular cloud. Because of the spherical focusing of the dust and gas inflows, the dust density 
increases to a level such that dust-dust collisions become important. Wh en dust grains col- 



l ide w ith e ach other they ma y agglo merate to form larger dust grains. iPrzygodda et al. 



(120031 ) and Ivan Boekel et al. I (120031 ) have repor ted d i rect o b servational evidence of grai n 
growth in circumstel l ar disks while, in addit ion, I Jural (119801). iMivake fc Nakagawal (119931 ) . 
Pollack et aFI Jl994h . IP'Alessio etaj] J200lh . and IPuUemond fc DominikI J2005h provided 
detailed calculations showing a strong tendency for dust grain growth when dust grains col- 
lide with each other. This agglomeration will increase the dust grain radius r^ while keeping 
the dust mass density pd constant. We will consider first how this agglomeration affects 
dust-gas collisions and then how it affects dust-dust collisions. 

Since the mean free path is much larger than the grain radius, the drag force on 
a dust grain due to colli sions with gas molecules is of the Epstein-type and given by 
(JLamers fc Cassinellilll999l ) 



drag 



-{ud - Ug)pgad\ d + {ud - u 



(12) 



where Cg is the gas thermal velocity, Ud is the dust grain velocity, a^ is the dust grain 
cross-sectional area, and Ug is the mean velocity of the gas (i.e., the fluid velocity). In the 
innermost Bondi region where flow is supersonic, we may approximate Cg ~ and work in a 
frame moving with Ug by defining Au^ = Ud — Ug. The dust equation of motion in this frame 
is thus 
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^d—7— = - {^Udf pgffd- (13) 

Defining ^ to be distance in the direction of dust motion so /S.Ud = d^/dt, Eq{T3] can be 
recast as 

^A«. = - {Audf aA . (14) 

d^ rud 

Integration gives 

AudiO = A«d(0) exp i-^/ldg) (15) 

wliere tfie dust-gas collision mean free path is 

kg = • (16) 

Pg(rd 

Since the dust cross-section and mass are given respectively by 

^d = Tirl (17) 

and 

m. = '-^ (18) 

where p^"* is the intrinsic density of a dust grain, the dust-gas collision mean free path can 
be expressed as 

A int 

ho = i^rd (19) 

which shows that dust agglomeration increases the dust-gas mean free path and so will tend 
to make dust collisionless with respect to gas. 

Let us now consider how agglomeration affects dust-dust collisions. We first note that 
the condition for dust-dust collisions to be significant is closely related to the condition 
for the dust to be optically thick: if / is the characteristic length of a configuration, the 
condition for collisions to be significant is pdCdl/rrid > 1 whereas the condition for the dust 
to be optically thick is QeffPd<^d/^^d > 1 where Qeff is an extinction efficiency parameter 
that depends on the ratio of the dust radius to the light wavelength. The dust-dust collision 
mean free path is thus 

,,, = J^ = 'J^U (20) 

PdC^d ipd 

so if, as argued in Paper I, the dust mass density pd has been enriched to be a significant 
fraction of the gas mass density pg, the dust-dust collision mean free path Idd will be the 
same order of magnitude as the dust-gas mean free path Idg- Agglomeration will thus tend 
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to increase both the dust-dust and dust-gas colhsion mean free paths, and furthermore will 
cause the dust to become optically thin. We will assume that dust grains agglomerate when 
the dust number density rid = Pd/T^d becomes sufficiently large for dust-dust collisions to 
occur and that this agglomeration results in an increase in r^ until the dust grains become 
collisionless and optically thin again. We will not attempt to follow the dynamics of the 
agglomeration process, relying instead on the analysis in the papers cited above. Our starting 
point then will be assuming the existence of collisionless dust grains exposed to star light, 
having radius r^ larger than in the ISM, and as discussed in Paper I, having a dust to gas 
mass density ratio substantially enriched compared to the 1% value in the ISM. 



4. Review: Neutral particle motion in a gravitational field 

For reference and in order to define terms to be used later in a more complex context, 
we first review the elementary problem of the motion of a neutral particle of mass m in the 
gravitational field of a star of mass M. The particle we have in mind could be a a dust grain 
with radius r^ sufficiently large to be collisionless over the distance from its starting point 
to the star. 

The equations governing the motion of this neutral particle are spherically symmetric 
whereas the motions of a charged particle in an azimuthally symmetric electromagnetic field 
are cylindrically symmetric. An axisymmetric magnetic field is assumed to exist in the lab 
frame and the z axis is defined by the direction of this magnetic field at the origin. Although 
the neutral particle trajectory is unaffected by this magnetic field, we nevertheless use the 
magnetic field coordinate system to define the lab frame. Depending on what is being 
emphasized, the lab frame will be characterized by either a cylindrical coordinate system 
{r, 0, 2;} or by a Cartesian coordinate system {x,y,z} so that x = rcoscp, and y = rsin0. 
Because the force is central, the neutral particle angular momentum vector L = mr x r is 
invariant and so the neutral particle moves in an orbital plane normal to L. The lab and 
orbital planes are sketched in FigjSJ The x axis of the lab frame is defined to be in the 
direction of the unit vector x = z x L/L and the y axis of the lab frame is defined to be in 
the direction of the unit vector y = z x (z x Li/L). The orbital plane is tilted with respect 
to the lab frame by an angle 6 about the x axis. The x' axis of the orbital frame is defined 
to be coincident with the x axis of the lab frame and the y' axis of the orbital plane is an 
uptilted version of the y axis of the lab frame. 

^ = corresponds to prograde motion in the lab frame (i.e., the neutral particle moves 
in the same sense as the toroidal current that produces the magnetic field B^ on the z axis), 
6 = 71 corresponds to retrograde motion in the lab frame, and 6 = tt/2 corresponds to a polar 
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Fig. 2. — Lab frame has Cartesian coordinates x, y, z and the magnetic field is axisymmetric 
with respect to the lab frame z axis. The orbital plane of a neutral particle is normal to the 
neutral particle angular momentum vector L which is tilted by an angle 6 with respect to 
the z axis. The orbital plane Cartesian coordinates are x',y' where the x' axis is coincident 
with the X axis. The neutral particle makes a circular Kepler, elliptical Kepler, or cometary 
orbit in its orbital plane (cometary orbit shown). 



cometary orbit 




y 
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orbit. For purposes of following the trajectory in the orbital plane it is convenient to use 
cylindrical coordinates p, rj defined in the orbital plane such that x' = pcosr] and y' = psinr/. 
The Hamiltonian for a neutral particle moving in its orbital plane can then be written as 

E = -mvl + ^ 21 

where 

L = mpVn, (22) 

the magnitude of the mechanical angular momentum vector, is an invariant positive scalar. 
The Kepler angular frequency at a reference radius a is defined as 



fio = VMG/a\ (23) 

The value of a is chosen to be the radius of the poloidal magnetic field axis. Normalized 
quantities are defined as 

p = p/a,T = not,Vp = ^ 

L = ^ H = ^., ^ ' 

Equation [21] can then be expressed in dimensionless form as 

H = S + 4-^. (26) 

2 2p2 p ^ ^ 

The last two terms depend on p and so constitute an effective potential 



1 
2p2 -p- 



Xip) = 77^ - -• (26) 



This effective potential depends parametrically on L which is a property of the particle 
and not the environment. Two different particles at the same position but having different 
values of L will have different effective potentials and so march to a "different drummer" . 
This "different drummer" concept will re-appear later in a more elaborate fashion when the 
motion of charged particles is considered. 

xir) attains its minimum value Xmm = —1/2L'^ at the normalized radius p = L^. A 
particle with energy equal to this minimum has Vp = and therefore has a circular orbit 
with angular frequency drj/dr = L/VLQima? = L/p^- Hence, if L = 1 the minimum-energy 
particle traces out a circular Kepler orbit with drj/dr = 1 and has an energy H = —1/2. A 
particle with energy —1/2 < H < cannot escape to infinity and so has a bounded elliptical 
Kepler orbit. The effective potential prescribed by Eq|26]for a particle with L = 1 is shown 
in Figint^a). 
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Fig. 3. — (a) Effective potential for a neutral particle having L = 1; (b) effective potential 
for a charged particle with appropriate values of canonical angular momentum and poloidal 
flux function. 



X 



2 
1 



r 



1 
2 








2 



3 



4 



X 



2 

1 



1 
2 








2 



3 



4 



-15- 



Reflection (pericenter) of a particle occurs when Vp = in which case Eqj25] gives 

p,e..cen*e. - ^^-^===. (27) 

Reflection at the pericenter can be considered to be the consequence of a potential barrier 
preventing the particle from accessing the region p < Ppericenter] the effective potential in 
the inaccessible region exceeds the total available energy. Thus an unbounded particle with 
L = 1 also has the effective potential shown Fig. [31(a), but unlike the bounded H = —1/2 
Kepler particle, the unbounded particle reflects from the pericenter potential barrier and so 
has a cometary orbit. 

In order for an incoming unbound particle to access a given p without being reflected at 
some larger radius, the condition that v"^ cannot be negative gives the constraint on angular 
momentum that 

l^ < 2p^H + 2p. (28) 

Since a particle with zero angular momentum will simply fall into the central object, in order 
for a particle to be both unbounded and able to access the radius p its angular momentum 
is constrained to lie in the range 

< L^ < 2p^H + 2p. (29) 



Solution of the equation of motion (JGoldsteinl Il950l ) shows that the orbit can be ex- 
pressed as 

1 _ 1- Vl + 2L2iJ-cos(?7-a) 

-p~ 12 ^^U) 

where a, which we call the clock angle in the orbital plane, is the angle between the symmetry 
line of the orbit (the line passing through the central object and the pericenter position) and 
the lab frame x axis (which is also the x' axis of the orbital plane). 

The Cartesian orbit coordinates x' = p cos r] and y' = p sin r] in the orbital plane (denoted 
by a prime to distinguish this plane from the lab frame) are 

™/ L^ cos rj 

1-a/i+2L2 Hcos{'q-a) /oi N 

\-^\+ll? Hcos{Ti-a) 

If the effective potential had a different shape, say the shape shown in Figl3]^b) with x ^ 
at large p, then a particle with H > could be trapped in one of the two minima of this 
effective potential. However, a particle coming from infinity would still be unbounded and 
would just reflect from some potential barrier. The inability of a static Hamiltonian system 
to trap a particle coming from infinity is independent of the shape of the Hamiltonian and 
results from the intrinsic time reversibility of Hamiltonian dynamics. 
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5. Comparison of gravitational/magnetic forces to Poynting-Robertson force 

and to radiation pressure 

The analysis in tliis paper is based on tlie assumption tliat tlie trajectory of cliarged 
dust grains results primarily from a competition between gravitational and magnetic forces 
with the possibility that in certain situations electrostatic forces and collisional drag can 
also be important. Two other types of forces, namely those due to the Poynting-Robertson 
effect and due to radiation pressure, also exist and so it is important to check to see if these 
additional forces need to be taken into account. This will be done by making a comparison 
with the nominal magnetic force on a charged dust grain. The magnetic force depends on the 
strength of the magnetic field, a quantity which has been estimated in self-consistent fashion 
in Paper III to be in the range 10~^ to 10~^ T (i.e., 0.1 to 10 mG) for a nominal YSO jet-disk 
system where the dust grains have coagulated to a nominal radius r^ =3 /im. This estimate 
of the magnetic field is in rough orde r of magnitude agreem ent with measurements rep orted 



by IChrysostomou et al. I (119941 ). by [Roberts et al. I (119971 ) and by lltoh et al. I (Il999l ). and 



also is in agreement with the expectation that the magnetic fields in a disk jet system should 
be much stronger than the nominal 10^^*^ T (i.e., 1 ^G) magnetic fields of the ISM. 

The radiation pressure acting on a dust grain at a distance r from a star with luminosity 
L is 

Prad = T ^Qrad{rd) (32) 

where Qrad{f^d) is the efficiency with which the photons are absorbed/reffected by the dust 
grain. This pressure results in a radial outwards force F^ad = PradO'd- If the dust grain radius 
is much larger than X^ad the wavelength of the radiation, then Qrad — 1 whereas if the dust 
grain radius is much smaller than the wavelength of the radiation then Qrad ~ {Xrad/rdY <^ 1- 
The nominal r^ ~ 3 /im dust grains assumed here are much larger than the nominal light 
wavelength and so Qrad ~ 1- 

Since the gravitational force Fg = mMG/r"^ is also in the radial direction, the force due 
to radiation pressure and gravity compete; the ratio of radiation pressure force to gravita- 
tional force on a dust grain is 

_ LQradO'd _ 3 L Qradi^d) /oo\ 

where EqsITT] and [18] have been used. Assuming r^ = 3 /xm, nominal luminosity L = Lq = 
4 X 10^^ watts, M = Mq, intrinsic dust density p**^* = 2 x 10^ kg m~^, and Qradijd) = 1 
gives a = 10~^ so radiation pressure can be ignored compared to gravitational force. 

The force on a dust grain due to the Poynting-Robertson effect is smaller by a factor v/c 
compared to the radiation pressure force, is in the toroidal direction, opposes the Keplerian 
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orbital motion vk = yMG/r, and so constitutes a drag force 



VK ^ QradTJ I MGL^ 

c Ac^ V r^ 



^Pn= Foa- = ^^i^\l^^. (34) 



The toroidal component of the magnetic force acting on a charged particle has magnitude 

Fmag = ZcVrB (35) 

where Z is the charge. The grains typically have non-circular trajectories with Vr being of 

the order of the Kepler velocity vk so F^nag ~ ZevxB. The ratio of Poynting-Robertson force 

to magnetic force is thus 

FpR Fg Vk ,„„. 

J—= J—^ — - (36) 

-^ mag ^ mag "^ 

Since the dust grains are assumed to be in a regime where they are acted on by magnetic 
forces which are at least comparable to gravitational forces, i.e., Fmag ^ Fq and since 
a ^ 1 and vk/c ^ 1 it is seen that the force due to Poynting-Robertson effect is negligible 
compared to magnetic forces and so the Poynting-Robertson effect, like radiation pressure, 
may be neglected. 



6. Electromagnetic particle Hamiltonian with gravity 

Hamilton-Lagrange methods are mathematically equivalent to the particle equation of 
motion and so describe all physically allowed orbits (e.g., cyclotron, drift, Speiser, etc.). 
Furthermore, because the Hamilton-Lagrange approach clarifies effects of spatial symmetries, 
deeper insight into orbital dynamics is obtained than provided by direct integration of the 
equation of motion. Direct integration nevertheless provides insight as well by providing an 
independent verification of the predictions of Hamilton-Lagrange methods. This two-pronged 
approach (Hamilton-Lagrange and direct orbit integration) provides a powerful method for 
examining particle motion in non-adiabatic situations. 

The Lagrangian of a particle with mass rrio- and charge q„ in the combination of an 
axisymmetric electromagnetic field and the spherically symmetric gravitational potential of 
a mass M central object is 

C = m^Ll + r'^^^ + vl 



+q„ (r^A^{r, z, t) + v,A,{r, z, t)) (37) 
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where V{r, z, t) is the electrostatic potential and a gauge with A^ = is assumed. The 
canonical angular momentum is 

dC 

P^ = —^ = m^r^(f) + q^rA^ (38) 

(90 

and, since Bpo; = V x (2nr)^^^(f) = V x f ^^0 j implies ip = 27rrA^, the canonical angular 
momentum can be expressed in terms of the poloidal magnetic flux as 

P^ = my^+^ij{r,z,t). (39) 

Zn 

Lagrange's equation P^ = dC/dcj) provides the important result that 

P^ = const., (40) 

i.e., P(^ is a constant of the motion because the system is axisymmetric. In the limit of a strong 
magnetic field, the second term in Eqj39] dominates the first and leads to the constraint that 
a particle orbit must stay very nearly on a surface of constant ip; this is the basis for particle 
confinement in axisymmetric toroidal fusion devices (tokamaks, reversed field pinches, and 
spheromaks). Any deviation of a particle from a constant ip surface is a consequence of finite 
rricr. When finite rUcr is taken into account, it is seen that the particle must stay within a 
poloidal Larmor radius of a constant ip surface, where poloidal Larmor radius means the 
cyclotron radius evaluated using the local poloidal field magnitude. Equation [32] niay be 
solved for to give 



The corresponding Hamiltonian is 



rrifj / 2 2i2 2\ . TA/ ,\ rricrMG 



H = ^ ivt + r'<l>' + vi]+q^V{r,z,t) - (42) 

By using Eqj4T]to substitute for in Eql42| the Hamiltonian can be expressed as 

H = n^{vl + vl) 

^(£V|g^ (43) 

We now consider situations where ip is time-independent and V = so the Hamiltonian 
reduces to 

2 y'r+v.)+ 2m.^2 v/HT^' ^ ^ 
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Since the Lagrangian does not explicitly depend on time, H = const, and the particle 
energy is conserved. In the q^tp = limit, P^f, reduces to the mechanical angular momentum 
P(j, = mr'^(p = mr x r-z = Lcos6 in which case the dynamics reduces to the neutral particle 
orbital mechanics reviewed in SecJH Thus, if q„ip = bounded orbits correspond to H < 
and circular Kepler orbits correspond to H having the minimum value of the effective 
potential well. Unbounded q^ip = orbits correspond to H >0. Note that p,/, = Lcos9 is a 
signed quantity, unlike L. 

If q^ip is finite then {P^ — qa'ilj{r,z)/27i) /2mo-r^ is the appropriate term which con- 
tributes to the effective potential. This term, called the Stormer potential, manifests a 
variety of qualitatively different spatial profiles depending on the relationship between P^ 
and q„ilj{r,z)/27T. These profiles are shown in FigJH for a sequence of decreasing values of 
P^. Very large positive P,/, gives prograde orbits similar to unmagnetized prograde cometary 
orbits and very large negative P^ gives retrograde orbits similar to unmagnetized retrograde 
cometary orbits; in both these cases the strong centrifugal repulsion at small r causes the 
particle to have an unbounded cometary orbit. 



Caption for Fig. H] 

Left: Plot of il^i^r, z) /ipQ v. r /a for 2; = with sequence of values of 27iP^/qipQ shown as 
dotted line. Right: Corresponding dependence of effective potential term (P^ — g-?/'(r, 2;) /27r) /r^ 
showing that potential wells develop at locations where 27iP^/qipo intersects ip{r, z)/ipo- These 
wells correspond to cyclotron motion if the intersection is away from the maximum of ijj and 
to Speiser orbits if the intersection is at or near the maximum of ip- A potential well at r = 
develops if P,^ = as seen in the sixth set of plots from top; this results in drain-hole orbits. 
Dotted line in right-hand second plot from top has the vertical scale multiplied by 100 to 
enable visualization of the outer minimum and the dotted line on the third plot from the top 
has the vertical scale multiplied by 2000 times. Note changes of scale in right-hand plots. 



-20- 



Fig. 4. — caption on previous page 
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On the other hand, if P^ and qa'ipir, z)/2it have comparable magnitude, complex effective 
potential structures can result. For example if at some location P^ equals Qa'ip^r, z)/2tt then 
[Pfj) — qfj'ilj{r,z)/27c) /2m^r'^ vanishes at this location, giving a localized minimum in the 
overall effective potential. 

If two separated positions exist where P<^ equals qai^ij", 0)/27r then two distinct minima 
exist, but if only one position exists where P^ equals go-'?/'(r, 0)/27r then only one minimum 
exists. The former situation occurs when P^ lies somewhere between and the maximum 
of ip and leads to cyclotron orbits with associated grad-P and curvature drifts; in this case 
EqJH] shows that the sign of (f) oscillates as the particle oscillates back and forth across the 
minimum of the effective potential, and the orbit is a cyclotron orbit. The situation of 
only one minimum occurs when the value of P^ is approximately the maximum of qrjip/2TT 
and gives Speiser orbits (0 has fixed sign and the orbit is paramagnetic). For a review of 
the distinction between the diamagnetism of cyclotron orbits (and associated grad-S and 
curvature drifts) and the paramagnetism of Speiser orbits see Appendix [Cl 

Yet another situation is where P^ = 0. Because t/^ ~ r^ at small r, this special 
case removes the singularity of (P,^ — go-'?/'(r, 0)/27r) /2m„r'^ at r = and so eliminates 
centrifugal force repulsion altogether. The P^ = case gives trajectories which spiral down 
towards the central object while crossing magnetic field lines; part of the magnetic force 
cancels the centrifugal force so all that is left is gravity and a residual inward magnetic force. 
This situation is completely different from either cyclotron orbits or Speiser orbits and has 
only been pr eviously discussed in the more limited context of two dimensional sit uations 



( jBellanl 120071 ). Finally, there is also the special situation discussed by ISchmidtl (119791 ) where 
Ptf, = —r'^qcr{27r)~^d{ip/r) /dr in which case the charged particle executes an axis-encircling 
cyclotron orbit. While possible in principle, axis-encircling cyclotron orbits will be not be 
considered here because they would correspond to particles having extreme energies (e.g., a 
cyclotron radius of many a.u.). 

Thus, there are five qualitatively distinct types of feasible trajectories depending on the 
relationship between P^ and qc,ip/2Ti. As labeled in FigJH and in order of descending signed 
value of the invariant P^ as shown by dashed horizontal line in left column of this figure, 
these are: 

1. prograde centrifugally dominated orbits [P^ much larger than the peak of qai'ir, z)/2tt)] 

2. Speiser orbits [P^ just grazes the peak of qa'>p{r, z)/2ti] 

3. cyclotron orbits [P^ well below the peak of qaip{r, z)/27t but much greater than zero] 
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4. P^ = orbits [which we will call "drain- hole" orbits for reasons to be discussed later], 
and 

5. retrograde centrifugally dominated orbits [P^ negative and much less than zero)]. 

In the above list, we have removed the constraint that z = so, for example, in case 
#3 (cyclotron orbits), the locations in the r-z plane where P<^ and q^tp{r, z)/27f are equal 
corresponds to a specific closed curve in the r-z plane; i.e., a specific ip iso-surface as shown 
in the projection of ip{r, z) at the top of Fig. [H 

The various possible values of P^ can be considered as the different "drummers" that 
dictate the effective potentials governing the motion of different particles located at the same 
position. A related e xample of this different "drummers" situation has been reported by 



Tripathi et al. I (120071 ) and involves two particles at the same location having velocities with 
equal magnitudes but opposite directions; the two particles have such extremely different ef- 
fective potentials that one particle is expelled from a magnetic flux tube (hill-shaped effective 
potential) whereas the other remains in the flux tube (valley-shaped effective potential). 



6.1. Mechanism for accretion of collisionless particles 

Accretion is the process of converting unbounded orbits (i.e., cometary orbits) into 
bounded orbits. Accretion of a collisionless neutral particle is clearly impossible if such a 
particle is governed by dynamics of a time-independent Lagrangian because converting an 
unbounded orbit into a bounded orbit would require changing the particle energy H and 
such a change is forbidden for a particle having a Lagrangian that does not explicitly depend 
on time. 

We now postulate an accretion mechanism as follows: photo- emission acts as an effective 
switch which alters the form of the Hamiltonian equation governing particle dynamics. The 
particle energy H and mechanical angular momentum mr'^cp do not change during the switch- 
ing, but after photo-emission has occurred, H and mr'^cf) become parameters in a different 
Hamiltonian system which has a different topography of potential barriers. For example, 
photo-emission can transform the neutral particle effective potential shown in Fig|31(a) into 
the charged particle effective potential shown in Figj3](b). 

The switching is postulated to occur when an incident neutral dust grain absorbs suf- 
ficient energetic photons from the star. The photon absorpti on causes the dust grain to 



photo -emit electrons and therefore become positively charged (JLed Il996l : ISickafoose et al. 



20001 ). The photo-emitted electrons become free electrons equal in number to the dust grain 



-23- 



charge Z . The photo-emission process has effectively caused the initial neutral dust grain 
to disintegrate into a single heavy positively charged fragment (the charged dust grain) and 
Z light negatively charged fragments (the photoelectrons). The motion of each fragment 
is governed by the Hamiltonian for a charged particle and this Hamiltonian is considerably 
different in form from the Hamiltonian that governed the neutral particle motion. 

The charge g^ of a dust grain charged by photo-emission is given by 

^^ - W^Hoton-W.^i (45) 



4vreor'd 
where Wphoton is the energy in eV of an incident photon that causes phot o-emission of a pri 



mary photo-electron and W^f is the work-function in eV of the material (jShukla fc Mamun 



20021 ). iLed (119961 ) has shown that the effective photon energy is Wphoton — 8 eV for nominal 
solar parameters and the effective work function of typical dust is W^j — 6 eV so that the 
energy of emitted photo-electrons is ~ 2 eV. 

Combination of Eqs. [18] and HS] show that the dust charge to mass ratio will be 

Qd _ Sgp (Wphoton — Wwf) . . „x 

which will be many of orders of magnitude smaller than the charge to mass ratios of electrons 
or ions. The number Z of charges on a dust grain will be 



ineord {Wphoton - W^f] 



(47) 



Charging a dust grain to qd takes a finite time interval, but for simplicity we will assume 
that this charging occurs at a single time defined as t = 0. Charging will not change either 
the instantaneous position or velocity of a particle. 

Photo-emission at t = therefore decomposes an incident neutral dust grain into positive 
and negative product particles each inheriting the same position and velocity at t = 0+ that 
the neutral dust grain had at t = 0_. Position and velocity can consequently be considered 
to be continuous functions at t = so the canonical angular momentum with which a newly 
formed charged particle is endowed is 

= Lf^cosO + ga'^(n,-2*)/27r 

where the subscript * denotes the value of a coordinate at the instant of charging, i.e., at 
t = 0. For simplicity we assume that photo-emission occurs when the distance between 
the incident neutral and the central object is at some critical spherical radius R^, so that 
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charging and the setting of t = occurs when the particle crosses the surface of the fictitious 
i?=K sphere, i.e., when r = r^, and z = z^, are such that 



2 I 2 

r^ +z. 



Rt 



(49) 



Depending on the value of L(j, the angle of inclination 6, the magnitude of g^, and the 
value of ip{r^, z^), all possible finite values of P^ can occur, including positive, negative, and 
zero. The r-z plane topography of the effective potential can change completely from what 
it was at t = 0_ because the centrifugal force potential p^/2mr^ responsible for the neutral 
particle potential barrier at small r is replaced by the Stormer term [P^ — q„ilj{r, z))^/27rmr^. 
Figure H] demonstrates that variation of particle mass and variation of the incoming orbit 
plane inclination angle 9 results in a range of p^, qd values and hence a range of P^ 
values corresponding to cyclotron, "drain- hole" , Speiser, or cometary orbits. If the new 
orbit is cyclotron, drain-hole, or Speiser, then photo-emission has prevented the particle 
from returning to infinity, i.e., the particle has accreted. Photo-emission changes the "rules 
of the game" by effectively erecting a new potential barrier which traps a previously unbound 
particle. The "old game" (i.e., neutral particle Keplerian motion as reviewed in SecH]) did 
not depend on particle mass or 9, but the "new game" does and leads to a mass- and 9- 
dependent sorting of incoming charged grains and their associated photo-emitted electrons 
into qualitatively different classes of orbits. 

This process whereby neutral particles enter a magnetic field from outside, become 
charged, and then become subject to magnetic forces is called 'neutral beam injection' in 
the context of tokamak physics and 'pickup' in the context of solar physics. Neutral beam 



injection is used r outinely for tokamak heating and current drive (jSimonen et al. 



1988 



Akers et al. 1 12002| ) . Pi ckup is important in the solar wind (JGloeckler fc GeissI l200ll ) , in 
plane tary atmospheres flHartle &: Killed l2006l ). and in producing source particles for comic 



rays ( Ellison et al. I Il998l ). However, to the best of the author's knowledge, charging of 



incoming neutral particles has not been previously proposed as a means for accreting matter 
around a star and instead accretion of matter around a star has always been argued to 
be the result of the viscosity of neutr al particles, i.e., collisio r is of n eutral particles with 
each o ther, as discus s ed fo r example in iLynden-Bell fc Pringld (Il974l ). IShakura fc Sunyaev 
(Il976l ). and iPringld (1198 ll ). The viscosity-based neutral particle accretion models suffer 
from not knowing what to do with the angular momentum of incident particles; this issue has 
motivated substantial work on developing the rather complicated non-linear turbulence-based 
magneto-rotational instability model as a means for transporting excess angular momentum 
outwards. In contrast, the model proposed here inherently accounts for angular momentum 
and so does not need any "add-on turbulence" to transport angular momentum outwards. 



Trapping via photo-emission has the remarkable feature that the special class of charged 
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particles created with zero canonical angular momentum will spiral all the way down to the 
central object. These P<^ = (drain-hole) particles falling towards r = are in what is 
effectively a loss cone in canonical angula r momentum space. The drain-hole particle mo- 



tion constitutes a gravity-driven dynamo (JBellanl 120071 ) because the accumulation of these 



particles near r = produces a radially outward electric field while their flow produces a ra- 
dially inward electric current (a dynamo is characterized by having opposed internal electric 
field and electric current). Since Jr = — (27rr)^ dl/dz (see Eqj5]), creation of this radially 
inward current which is symmetric with respect to z implies creation of an ant i- symmetric 
function /(r, z) which in turn implies creation of an ant i- symmetric toroidal field B^ (see 
EqJl]). Equally remarkable, particles for which P^ is near the maximum of q„tp/27r develop 
Speiser-type paramagnetic orbits in the vicinity of r = a, 2; = and so can constitute the 
toroidal current that produces the poloidal flux (see discussion of Speiser orbit paramag- 
netism in Appendix [n]). The creation of Speiser-orbit particles is concept ually similar to 



toroidal current drive in a tokamak via tangential neutral beam injection (jSimonen et al. 



19881 ). Because the drain-hole and Speiser dynamos ar e both axisyni metric, both violate 



the essential claim of Cowling's anti-dynamo theorem (j Cowling! Il934l ) that axisymmetric 
dynamos cannot exist. This violation is not a problem because Cowling's theorem is based 
on MHD and so does not take into account drain-hole or Speiser orbits. 

The Hamiltonian for an incoming neutral dust grain of mass m„ can be written as 

_ rrinvl rrinvl mnvl ninMG 

This neutral dust grain absorbs energetic photons at t = 0, photo-emits Z free electrons, 
and consequently becomes positively charged with a charge of Z. The mass of the neutral is 
related to the mass ?7i_|_ of the positively charged dust grain by ?n„ = ttlj^ + Zm,, where rrif. 
is the electron mass. Prior to this charging process, EqlSO] can be written as 

TT ^ {m++Zme:)vl . {m++Zm.e)vl 

, {m++Zme)vl _ {m++Zme)MG V '' 

+■ 2 V^F^T^ 

so the positively charged dust grain with mass m+ and the Z electrons can each be thought 
of as executing identical neutral-type orbits before photo-emission occurs. 

At the instant before charging, the neutral particle mechanical angular momentum is 

p^ = mnrlcf)^. (52) 

At the instant after charging the newly created positively charged dust grain and its asso- 
ciated photo-emitted electrons all have the same values of r^ and 0* . The canonical angular 



-26- 

momentum of the positively charged dust grain will therefore be 

P^ = m^rlcj)^ + Zeip{r^, z^)/2'n (53) 

and the canonical momentum of each associated electron will be 

PI = rrierl^^ - e^{n, z^)/2tt. (54) 

The initial neutral dust grain will be called the "parent" particle while the positively charged 
dust grain resulting from photo-emission and its associated Z photo-emitted electrons will 
be called "child particles" that are "siblings" of each other. The set of child particles re- 
sulting from the charging of a specific neutral dust grain will be called a "family". The 
canonical momenta Pt and Pf are now the appropriate orbit invariants for t > whereas 
the mechanical angular momenta m„r'^<p of the individual siblings will not be invariant for 
t > 0. Although the mechanical angular momentum of an individual sibling is not conserved, 
the total mechanical angular momentum of the family is conserved since summing Eqj53] and 
Z times Eqj5l gives [P4,]f^jj,iiy = -^/ + ^^4> where [p<i,]f^^iiy = mnrlcj) is the sum of the me- 
chanical angular momentum of the charged dust grain and all its sibling electrons. Because 
the siblings can physically separate from each other, the mechanical angular momentum 
[p^ ,^^i is not a locally defined quantity and so is not a constant of the motion of either 
a single particle or, as in ideal hydrodynamics, of a fluid element. However, the mechanical 
angular momentum of the entire system is conserved because the angular momentum of each 
family is globally conserved. 

The kinetic energy of each sibling at the instant before photo-emission is the same as 
the value at the instant after photo-emission. If H is decomposed into the contributions from 
the various siblings, it is seen that each sibling's H^^ is the same before and after photo- 
emission. Assuming zero electrostatic potential for now, but allowing the child particle to 
be at arbitrary z, the Hamiltonian of each sibling is 






rriaVt 



2 ' 2 

r 

2mcrr'^ 



rrigMG 



Using EqlHlwe now define 



{ujccj) = (56) 

rrirr 



as the spatially-averaged cyclotron frequency over the area bounded by the poloidal field 
magnetic axis. We now normalize all quantities to appropriate combinations of a and Qq, 
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the Kepler angular frequency at a prescribed by Eqj23l The normalized time, cylindrical 
coordinates, velocities, and magnetic flux are thus 



in which case EqJSSl becomes 



T 


= 


rtot 


f 
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r/a 
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z/a 


Vr 
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Vr/aVt^ 


Vz 
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Vz/aVlQ 
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L^/m^a^VtQ 
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f'^d(f)/dT = L cos 6 

%l){r,z) _ ip{r,z) 
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mQ.QO^ 



(57) 



H 



"r J ^ 

2 ~'~ 2 



2^2 



\/r'^+z 



The effective potential is now 



X{.r,z) 



mechanical 
Lcos6' 



magnetic 



(5^ 



-^^ [^(r„2;,)-7/'(r,^)] 



V 



(59) 



2 f2 



Vr^ + z^ 
gravitational 
where the mechanical, magnetic, and gravitational contributions have been labeled. Before 



photo-emission, the mechanical angular momentum is invariant so 



I ' * , Z-j. 



Lcos9 is just 



the normalized mechanical angular momentum that the parent had when it was at infinity. 
Invoking Eqj29t it is seen that 



0<L< X 2R'iH + 2R^ 



(60) 

is required since incident neutral dust grains with mechanical angular momentum outside 
this range would have reflected at larger radii than R^, and so would not have been able to 
access the radius R^.. 
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The normalized canonical angular momentum with which a typical sibling charged par- 
ticle is endowed is 

p^ = Lcose + ^^{n,z,). (61) 

The parameters underlying Figj3]^b) can now be understood. This figure is a plot of x(^) z) 
V. f ioT z = where P^ is calculated for the situation where (uca) /^o = 40, n = 0.5, ^* = 0, 
6 = 0, and L = 1. Charging of an L = 1 dust grain via photo-emission causes the effective 
potential governing the dust grain to change from the form given in Fig|3]^a) to the form 
given in Fig|3](b). 

When the magnetic term in Eg J59] becomes comparable to the mechanical term or much 
larger, orbital dynamics for the siblings become very different from the orbital dynamics 
of the neutral parent that existed before photo-emission. Various orbits can occur for the 
siblings. Because of the complexity of these three dimensional orbits, we will first consider 
orbits confined to the z = plane and then generalize to fully three dimensional orbits 
ranging over finite z. 



6.2. Distribution of Cometary, Speiser, Cyclotron, and Drain-Hole Orbits 

As reviewed in SecJH neutral particles orbits are degenerate with respect to their orbital 

plane inclination angle 6 (see Figj2]). However, once particles become charged, they are no 

longer restricted to an orbital plane, and furthermore, as seen from Eqj59l the effective 

potential of a charged particle has a strong dependence on the value of 9 that its parent 

particle had. This dependence was manifested in the discussion of Fig. H] where it was 

noted that particles with \2'7iPif,/qa-i'o\ ^ 1 are essentially unmagnetized and have Keplerian 

cometary orbits, particles with 2'irP^6/qrj'ipo — 1 have Speiser orbits, particles with ^ 

27rP(^/go-^o ^ 1 have cyclotron orbits, and particles with 2nPtf)/qaipo ~ have drain-hole 

orbits. This discussion can be made more quantitative by defining A = 27iP^/q^ipQ] note 

that A corresponds to the horizontal dashed lines in the left hand column of FigJH Equation 

[6T] can then be recast as 

— 2Q 

A=—^Lcose + ^in,z,). (62) 

(^ca) 

Thus particles with |A| ^ 1 have Keplerian cometary orbits, particles with A ~ 1 have 
Speiser orbits, particles with ^ A ^ 1 have cyclotron orbits, and particles with A ~ 
have drain-hole orbits. 

Assuming H <^ I and ^=k ~ 1, Eql60] implies that only particles with < Z < ^/2 
can access a given location. Because there will be a distribution of all possible L's within 
this allowed range, we consider a particle with the mean of these allowed values as being 
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representative and so assume that L = v^/2 is the normahzed angular momentum of this 
representative nominal particle. 



Since (oOca) = Qa {Bz) /Tn„, {Bz) = ipo/'Ka'^, and fio = \/ MG jo? this nominal particle 
will have 

k = Kco^e + i){r,,z,) (63) 

where 

K = 64 

parameterizes the competition between gravitational and magnetic forces. Using Eq. H6]to 
give the charge to mass ratio it is seen that 



3£:o {Wphoton - W^f) tpo '^' 

thus K increases when r^ increases as a result of dust grain coagulation. 

Speiser and drain hole particles occur when gravitational and magnetic forces are com- 
parable in magnitude, i.e., when K is of order unity. For a given star mass M, poloidal 
flux magnetic axis radius a, and magnetic flux ijjo, this means that Speiser and drain hole 
particles will occur when coagulation has caused the dust grains to have a certain critical 
radius which is of order 



' (2aMG)'/'V ""Pd 

If Vd ^ rj'^^ then gravity will dominate and the dust grains will behave like neutral particles 
whereas if r^ -C r^**** then magnetic forces will dominate and charged dust grains will mainly 
have cyclotron orbits. Since coagulation causes r^ to increase monotonically, there should 
always be some time when r^ ~ r^*"** and K is of order unity. This argument indicates that 
the dust-driven dynamo mechanism should take place as a well-defined temporal stage in 
the accretion process; before this stage r^ is too small and after this stage r^ is too large. 

Since il){f^:,z^..) ranges between and 1, let us consider the nominal situation where 
V'(^*, ^*) = 1/2 in which case 

Anon, = K cose + -• (67) 

If i^ ~ 1/2, Speiser particles result for cos 6' = 1 (i.e., neutral parent was prograde) and 
drain-hole particles result for cos6 = —1 (i.e., neutral parent was retrograde ). li K <^ 1/2, 
the orbits will be cyclotron. Finally ii K ^ 1/2, the orbits will be cometary if cos^ is not 
close to zero. The categorization implied by EqlHTlis shown schematically in Fig. O 
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Fig. 5. — Distribution of orbits as function of K, 9. Radius K is proportional to r^. Prograde 
orbits have ^ = 0, retrograde orbits have 6 = n, and polar orbits have \6\ = 7r/2. 
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6.3. Light-weight particles {K ^ 1) 

6.3.1. Generic accretion mechanism 

Potential barriers occur at locations where x(^? ^) > H. We first consider motion of a 
sibling particle constrained to stay in the z = plane (i.e., the particle starts with Vz = 
and no forces exist that push it off the z = plane). In this case the effective potential is a 
function of f only and is 



/ 



mechanical 
Lcos^ 



magnetic 



\ 



+ 






[^(f„0)-^(f,0)] 



X{r,0) 



\ 



J 



2f 






gravitational 
f 



(68) 



If x(f, 0) exceeds H at some radius f > f^,, the particle becomes trapped within a finite 
extent region as indicated in Fig|3]^b) or equivalently by the third and fourth rows, right 
hand column of Fig. |H If f^, is small or large compared to unity so r^, is not near the peak 
of ip, then ipi^f^, 0) ^ 1. Because the particle is assumed to be light-weight (i.e., r^ -C r^"*), 
its average-field cyclotron frequency (ujca) will be much larger than the Kepler frequency Qq. 
Since L is of order unity, the light-weight particle will have \{ujca) /2fio| ^ -^ in which case 



max 



LcosO 



2Qo 



n, 0) - 7/;(f, 0)] 



2f 






max 



2f2 V 2^0 



V^(r,0) 



so the peak of x(^?0) will occur where '(/'(r,0) takes on its maximum value, namely unity. 
The maximum of x{^i 0) fo^' ^ light-weight particle is thus 



y(l ^)-- ^i^ 



(69) 



where the —1 term comes from the gravitational potential at r = 1, 2; = 0. This gives the 
necessary condition for trapping light-weight charged dust grains to be 



\{^ca)\ 

fin 



>2V2. 



(70) 



Because a typical light-weight grain has Kli^co-)! ^ 2v^fio, a light-weight grain confined to 
the 2; = plane will become trapped in a finite-sized region upon being charged, i.e., it will 
have accreted. The same will be true for the associated sibling photo-electrons since they 
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Fig. 6. — An incident neutral particle {H = 0, Ppericenter = 0.1, ^ = 30, a = 0°) becomes 
charged due to photoemission of electrons at i?* = 2. The child particle mass is such that 
(ujca) /^o = 200 and the child particle becomes magnetically trapped, staying within a 
poloidal Larmor orbit of a constant ip surface. In this example, the child particle is mirror 
trapped and so cannot enter the strong magnetic field region at small f. (a) x — y plane, 
charging occurs where orbit abruptly changes, poloidal field magnetic axis shown as dashed 
circle (b) f — z plane showing magnetic mirroring of child particle at large magnetic field 
(poloidal flux contours ip{f, z) shown as dashed lines). The orbit the parent neutral particle 
would have continued to have if it had not become charged is shown by dotted line (both 
projections). 
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also have K^Jce)! ^ 2v^f2o. Figure [6] shows a direct numerical integration of the equation 
of motion demonstrating this basic accretion mechanism in three dimensions: a hght- weight 
neutral dust grain disintegrates at a certain location into an K^iJco-)! ^ 2v^fio positively 
charged dust grain (there would also be Z associated photo-electrons which for clarity are 
not shown in the figure but would also have cyclotron-type orbits). An actual dust grain 
would start with an infinitesimal energy < ff -C 1; the calculation here uses /f = as 
representative of this infinitesimal H since the difference between an orbit with H = and an 
orbit with infinitesimal H is insignificant at any finite distance. The newly created charged 
particles are trapped to the vicinity of a ip{f, z) = const, poloidal flux surface (poloidal flux 
surfaces are shown by dashed lines in FiglHt^b)). Because of fi conservation, the charged 
particles can also be mirror-trapped, so while on the the constant ip surface, they reflect 
from regions of this surface where the magnetic field is strong. Figure Wl^a) plots the time 
dependence of p^ for the particle shown in Figj6l 

Figure[71^b) plots the canonical angular momentum P^ (solid line) and the kinetic/potential 
energies (dashed lines labeled 'KE' and 'PE' in figure). It is seen that p,^ is conserved before 
charging whereas P^ is the conserved quantity after charging. Also, the total energy (kinetic 
-|- potential, dashed line labeled 'Tot' in figure) remains zero. Strictly speaking, this plot 
should be considered as referring to the neutral dust grain until charging, and then to the 
charged dust grain after charging so the jump in P^ at the charging time seen in the figure 
does not violate the requirement that P^ is a constant of the motion for a specific particle. 

Figure M shows the three dimensional orbit of a light-weight charged dust grain with 
slightly different parameters so that it is not mirror trapped. The derivation of the non- 
dimensional equation of motion used here is given in Appendix [Dl 



6.3.2. Width of light-weight particle trapping well and relation to cyclotron orbits 

If KcUco-)! /2f2o ^ \p^{f^:,z^:)\, the magnetic term in EqJ59] dominates the mechanical 
term as soon as f deviates slightly from r^,. This implies existence of a narrow trench-like 
potential well with minimum very close to r^,. The effective potential shown Figl3]^b) has 
such a trench; this situation involves a particle confined to the z = plane and the trench 
is at f = p = 0.55. This situation is also evident in the third and fourth rows, right hand 
column of Fig. HI If | (ujca) \ ^ ^o fhe gravitational term is completely overwhelmed by the 
magnetic term so the trench bottom in the z = plane is where 

V^(^) = r^P</>(n,^,) +V^(r,). (71) 
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Fig. 7. — (a) Mechanical angular momentum p^ v. time r and (b) kinetic energy (KE), 
potential energy (PE) and canonical angular momentum P^ v. time for the calculation 
shown in Fig. [6l Mechanical angular momentum p^ is conserved before charging, but 
oscillates after charging; canonical angular momentum P^ is much larger than mechanical 
angular momentum because of strong magnetic field and is conserved after charging. 
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Fig. 8. — Same parameters as Fig. El except -R* = 0.8. Charged particle is now not mirror- 
trapped. 
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Taylor expansion of iplr) near r* gives 



d± 

(72) 



+I('=-'=.)HS --+■■■ 

If f^, is not close to unity, then ^(f^,) is not close to its maximum value so the leading term 
in the Taylor expansion is the one involving dip /dr. Using EqI72]to substitute for ■il>[f) in 
EqJTT] gives the trench bottom to be at 

so the trench bottom is close to f* because K^iJccr)! ^ fio is being assumed. If f* < 1 
then ipil) / df) ._._ is positive and vice versa since if) has its maximum value at f = 1. The 
trench bottom will thus be outside of f* if f^, < 1 so {ujca) [dip/df) ___ is positive and vice 
versa if f,,, > 1. The sign of ijj{f) — ^(f*) oscillates as the particle bounces back and 
forth in the trench. Using EqE] expressed in normalized variables, and noting that P^ = 
p(j,{r*, ^*) + 'ipir*, ^*) {^ca) /2fio it is seen that the azimuthal velocity 

2% 

dr ^^ 



"^ -W^(^") 



P0(f*X) -^[m-^{r,)] ^ '* 



has an oscillating polarity. The combined oscillation of r and dcp/dr c orresponds to the 



particle tracing out Larmor orbits with gyro-center at the trench bottom (jSchmidtlll979l ). 



To summarize: In the \{uca)\/^o ^ 1 situation (i.e., light-weight particles) Eq. [59] 
provides an effective potential whereby the magnetized charged particle is c onfined to the 



vicin ity of a constant ip surface, just like a charged particle in a tokamak (JRome fc Peng 



19791 ). The particle motion over the constant i/j surface can be understood as a sum of 
parallel to B motion, cyclotron motion, and particle drifts (curvature, grad B, etc.) as given 
by Eq JC5[ Furthermore, regions where fiB is large can constitute an additional potential 
barrier (i.e., magnetic mirror) that prevents those subsets of particles having inadequate 
velocity parallel to B from accessing the entire constant ip surface. Because of magnetic 
mirroring by ^VB forces, particles in these subsets are confined to the weaker magnetic field 
regions of a constant ip surface as seen in Figj6](b). 



6.4. Speiser orbit particles (KcosO ~ 1/2) 

When f^, ~ 1 and ^^, ~ 0, photoemission occurs near the peak of iplf, z), i.e., where V^ 
0; see Fig. H] second row from top where P^ is just grazing the peak of ip. The linear term 
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in the Taylor expansion in Eql72] is therefore neghgible. Since d'^ip/dr'^ is negative near the 
maximum of if), EqJ72] becomes 



m^i'ir.)- \{f-f.f 






ion |74| then reduces to 


^^ P0(--*^O) +i^ {f-f.f 





(75) 



dr 






(76) 



and for p^if^., 0) being positive (i.e., parent particle was prograde), dcp/dr is always positive. 
This corresponds to Speiser-type orbits because when the particles bounce back and forth 
across the peak of ip, they are bouncing back and forth between regions where the poloidal 
magnetic field ~ dip/dr changes sign. As discussed in Section [C73] of Appendix [C | this results 
in paramagnetism, i.e., positively charged particles moving in the positive direction and so 
producing rather than opposing a B^ field. Creation of Speiser-or biting particles sustains the 
poloidal magnetic field against losses and will amplify an initial seed poloidal field; creation 
of Speiser particles therefore constitutes a dynamo for driving toroidal current. 

Figure [9] shows the creation of a Speiser orbit by photo-emission charging of a neutral 
particle near the poloidal field magnetic axis. Figure [9](a) shows that the orbit is paramag- 
netic (i.e., particle moves in positive direction) while Fig. [9](b) shows that the orbit involves 
repeated reflection from the interior of a poloidal flux surface in the manner discussed in 
Section IC.3I of Appendix O 



6.5. Drain-hole particles (i^cos^ ~ —1/2) 



Charged particles born with P^ = are called 'drain-hole' particles because they behave 
as if they are going down a drain. The propertie s of drain-hole particles restricted to the 
z = plane were briefly examined in iBellanI (120071 ): here the more general 3D situation will 



be considered. Using Eqj6T] and 



I ' *i z^ J 



LcosO 



-- L cos 9 the P,^ 



condition corresponds to 



2fin 



i^ir^ 



(77) 



implying that cos 6 is negative in which case the parent particle must have been retrograde. 
The effective potential (see Eq l59|) for the P^ = class of particles reduces to 
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Fig. 9. — Speiser orbit resulting from parent with H = 0, ppericenter = 0.95, 6 = 18°, a = 0°. 
Charging occurs at -R* = 1.2 and the child particle is a positive particle with ujca/^o = 10; 
(a) the x-y plane orbit is counter-clockwise corresponding to paramagnetic motion; (b) the 
f-z plane orbit involves the particle continuously reflecting from the interior of a toroidal 
flux tube. 
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which has a funnel (i.e., drain-like) profile near f = 0, z = due to the second (gravitational) 
term and a hill on the funnel side wall with peak near f = a, z = due to the first (Stormer) 
term. A particle initially on the hill (i.e., near the poloidal field magnetic axis) will fall down 
the hill into the drain-like funnel; see right hand column, second row from bottom in Fig. H] 
for plot of first term in EqJTSl Thus, no matter where a P^ = particle starts in f, z space, 
it eventually follows a spiral path down to f = 0, z = 0; no centrifugal force will ever push 
it back outwards because the first term in Eq|78] has no singularity at f = (recall that 
?/' ~ f^ for small f, z). The sense of this downward spiraling trajectory will be in the — (tUco-) 
direction as shown by Eql74l Since t/^ ~ f^ for small f, z it is seen from Eq|74]that drain-hole 
particles have a limiting angular velocity 

hm ^ = - %1 hm (^) = const. (79) 

Combination of Eqs. [29] and [771 show that drain- hole particles can only be created if the 
accessibility condition 

^^ (^(n, z,)Y < {2p'H + 2p) cos^ 9 (80) 

is satisfied, a condition that {uca) /^o not be too large. Since /f ~ is assumed, cos6' ~ —1 
for drain-hole particles, and since p = \/r^+z^ is required to be larger than the pericenter, 
this condition becomes 

\y\'*'>'^*)\ ^ "/ \ ~ y ^ Ppericenter ■ \^^) 



Because EqJHT] requires (ujca) /^o to be small, drain-hole particles, like Speiser particles, 
result from dust grains with grain radii consistent in order of magnitude with EqJMl 

Particles with P^ exactly zero ( "perfect" drain-hole particles) fall down the gravitational 
potential all the way to the central object at the origin f = 0, z = 0. Particles that are not 
quite perfect drain-hole particles will have small, but finite P^ and so will reflect when close 
to the central object. 

Figure [To] shows a numerical calculation of a drain- hole particle orbit with H = 0. The 
solid line in Fig. [TO](a) shows the projection of the drain-hole orbit in the x-y plane. The 
orbit the neutral particle would have had if it had not encountered any photons and so 
remained neutral is shown as a dotted line. Figure [TUl^ b) shows the projection in the f-z 
plane with surfaces of constant ip indicated (the dotted line again shows the orbit the neutral 
particle would have had if it had not encountered any photons). The drain-hole particle has 
a retrograde orbit (clockwise sense resulting from its angle of inclination 6 > 90°). Figure 
[TlTa) shows that the mechanical angular momentum is not constant after charging while Fig. 
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Fig. 10. — Drain-hole particle {H = 0, Ppericenter = 0.2, 9 = 170°, a = 0, .R,^ = 0.8) falls 
across magnetic field all the way to the central object. This is a heavy particle (dust grain) 
and has Uca/^o = 1-6; (a) shows orbit projection in x-y plane, (b) shows projection in f-z 
plane. Dotted line shows trajectory parent would have continued to have if it had not become 
charged. 
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Fig. 11. — (a) Mechanical angular momentum p^ for drain- hole particle is not conserved when 
particle becomes charged, (b) canonical angular momentum P^, kinetic energy (KE) and 
potential energy (PE). The magnitudes of the kinetic and potential energy increase without 
bound as the particle falls towards the central object. The canonical angular momentum is 
conserved and is near zero upon charging. 
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[TT](b) shows that the canonical angular momentum remains constant at zero after charging. 
Figure [TTjfb) also shows how the magnitudes of the potential and kinetic energies increase 
without bound as the particle descends towards the central object while the total energy 
stays zero (kinetic, potential and total energies shown as dashed lines). 

When drain-hole particles approach the central object, the gravitational term in EqlTHl 
dominates (recall that ■?/; ~ f^ at small f and near z = 0). It therefore makes sense to use 
spherical coordinates in this region in which case the Hamiltonian is approximately 

1 /dRV 1 



~ - ^ ^ - - (82) 

2\dT ) R ^ ' 

where R is the spherical radius and ff ~ has been assumed. Equation [H2] shows that the 
free-fall velocity scales as 



dR 



dr 



w f*^^) 



and particle flux conservation over a spherical surface AnR"^ shows that ATTR?ndh{R)dR/dT = 
const, where ndh{R) is the density of drain-hole particles. Thus, if the incoming drain- 
hole particles do not accumulate, spherical focusing combined with the accelerating free-fall 
velocity shows that the density of drain-hole particles scales as 

""''^"^^ ^ WMIdTr ^ Wr-- ^''^ 

Accumulation of the drain-hole particles in the vicinity of the central object will also 
increase the density of drain-hole particles with time. There is thus both a temporal increase 
and a geometrically-induced increase of the drain-hole particle density as R decreases. Since 
the sibling electrons were left stranded at large f, what results is the establishment of 
a large positive charge density near the central object and an equal-magnitude negative 
charge density at large f. The flow pattern of the drain hole particles and the location of 
the stranded electrons is sketched in FigJT^ Eventually the positive space charge near the 
central object becomes so large that it produces a repulsive electrostatic electric field that 
balances the gravitational force acting on any additional drain hole particles. The large 
positive potential near f = will tend to drive axial electric currents flowing away from the 
z = plane resulting in the loss or neutralization of some of the drain-hole particles. The 
axial electric current could result from attraction of electrons near f = towards the 2 = 
plane or from expulsion of positive particles away from the ^ = plane. Either electron 
attraction or positive particle repulsion will deplete the positive space charge density near 
f = 0. There will then have to be a replenishing flow of additional drain hole particles into the 
f = region to compensate for this depletion of positive space charge. Being very low mass. 
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the stranded electrons have very small Larmor orbit radius and so are constrained to stay 
essentially right on the poloidal flux surface on which they were photo-emitted (see Sees J6. 31 
and l6.3.21) . The electron flow is thus at a much larger \z\ than the drain-hole particle radially 
inward flow which is concentrated near the 2 = plane. The vertical separation between 
the respective radially inward flows of positive and negative particles means that bipolar 
toroidal magnetic flelds will be generated in the interstitial regions between the electron flow 
and the drain-hole flow (see positive and negative B^ regions in Fig IT2l) . 

If no electric current is allowed to flow, the situation is like a free-standing battery not 
connected to any load, i.e., a situation where there is a voltage differential across the battery 
terminals, but no current flows. However, if bipolar axial currents are allowed to flow, then 
the situation is like a battery connected to a load and the resulting radially inward drain- 
hole particle current in the ^ = plane is like the internal current in a battery. The overall 
current flow pattern sketched in FigJT2l results from a combination of drain hole particle and 
electron motion. This pattern is sketched in Fig{T3]as a conventional electric current. The 
geometry of the current flow pattern and electromotive force driving this current is identical 
to the geometry and flow patterns in the l abora tory conflguration s imulating astrophysical 



jets described in iHsu fc BellanI ( 12002I . 120051 ) and lBellan et al.\ (120051 ). The electric fleld due 
to the drain hole particles corresponds to the electric fleld produced by the capacitor bank 
used in the laboratory expe riment. This geometry and symmetry is also identical to that 



proposed by iLovelacd (119761 ) , the only difference being the means by which the radial electric 
fleld is produced. The magnetic flelds in the lab and astrophysical plasmas have the same 
toroidal/poloidal topology. 

The drain-hole current is thus powered by gravity and has Jr radially inward with 
E^ radially outward so that J ■ E is negative, consistent with the condition for a dynamo. 
The drain-hole particles have retrograde motion so their mechanical angular momentum is 
negative. This negative mechanical angular momentum is removed by the braking torque 
r X F = {rf) X [J^r x Bzz) = —rJrBzZ which is positive since Jr is negative and Bz is 
positive. 

The radially inward current is symmetric with respect to z. This property provides 
enough information to determine the symmetry properties of I{r,z). The radially inward 
drain-hole current means that J^ < and J^ = in the z = plane. Since Eql5] shows 
that Jr = —{27Tr)~^dI/dz and Jz = {27fr)^^dl/dr the condition Jz = means I{r,z) 
must vanish in the z = plane. Furthermore / must be an odd function of z in order for 
Jr = —{2TTr)~^dI/dz to be flnite in the 2; = plane. Finally dl jdz should be positive in 
order to have Jr < 0. Thus, /(r, z) should be positive for z > and negative for 2; < so that, 
as sketched in FigJT3| there will be a bipolar axial current flowing along the z axis outwards 
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Fig. 12. — Drain hole dust grains fall across poloidal field lines towards central object leaving 

behind stranded electrons which are confined to poloidal flux surface on which they are born. 

Drain hole particles accumulate near central object creating large positive charge there. This 

repels positive particles (drain hole particles, ions) to flow axially away from z = plane 

and also attracts stranded electrons which can flow on poloidal flux surface. The result is a 

clockwise poloidal current flow pattern in upper-half r-z plane, giving a positive B^p in region 

linked by poloidal current and a negative B^ in lower-half r-z plane where poloidal current 

flow is counter-clockwise. 
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Fig. 13. — Flow of conventional electric current for drain hole particles and their associated 
stranded electrons. The electric field on the z = plane is radially outwards while the 
current flow is radially inwards so J ■ E is negative, indicating that the infall of the drain- 
hole particles constitutes a dynamo. The J x B force (which is essentially due to the gradient 
of -B? and which is strongest at small r) drives a bipolar axial jet. 
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from the z = plane. The accumulation of drain-hole particles constitutes the engine 
that drives the poloidal electric current that drives the astrophysical j et. The ^-symmetry 



of ijj(r , z) and z-antisymmetry of /(r, z) has been noted previously by iFerreira fc Pelletier 



(119951) 



Electromagnetic power flow from this dynamo can also be interpreted in terms of the 
Poynting flux S = E x B/yUo- Azimuthal symmetry applied to Faraday's law shows that E^ is 
zero for a steady-state situation in which case the 2;-component of the Poynting flux reduces 
to 5*2 = ErB^I iiQ. Because Er and B^ are both positive for 2; > whereas E^ is positive 
while B^ is negative for z < 0, it is seen that Sz is positive for 2; > and negative for 2; < 0. 
Thus, the Poynting flux associated with this dust-driven dynamo injects energy into bipolar 
astrophysical jets flowing normally outward from the z = plane. 

One can ask just how close to exactly zero P^ has to be in order for a particle to behave 
as a drain-hole particle. Exact P,/, = would enable a particle to spiral down all the way to 
the center of the central object, an obviously unrealistic situation because the particle would 
vaporize as it approached the stellar surface. A more realistic question then is how small 
does P^ have to be in order for a drain-hole particle to fall to some specified normalized 
radius Rsmaii that is much less than unity. Rsmaii would presumably be of the order of the 
radius at which the astrophysical jet starts and so would be of the order of the thickness of 
the accretion disk or somewhat smaller. Since the dimensionless form of EqHH is 

H= ^(t^.'+t^D +Xir,z) (85) 

where the effective potential is 

Uf,z) = ^(p. - i^^(f,z)^ - , ^ (86) 

and since -ff ~ 0, the turning point for a drain-hole particle is where x(^? z) ~ 0. Because 
ijj{f, z) ^ at small f, the inner turning point will therefore be where P| = 2r^ j\/f'^ + zP-. 
Assuming that the inner turning point is at f ~ Rsmaii and 5 ~ 0, the inner turning point 
is where -P| = 2Rsmaii- A sufficient condition for assuming P^ ~ is thus P^ < 2Rsmaa, i-e.. 



- \/2Rsmall < LcOSe+ — ^ ^{r^, Z^) < \^2Rsmall (87) 

and particles satisfying this condition will fall to a normalized radius R < Rsmaii- For given 
L, {oJca) /2VLq, and Tp{r^, z^) this corresponds to a narrow range in 6 centered about the angle 
at which P^ equals zero exactly. Equation EZI can be expressed as cos {6 + A9/2) < cos 9 < 
cos [6 - A6'/2) where 

Zcos(^ + A^/2) = -y/2R^i-^ij{f,,z,) 
Lcos{e-Ae/2) = ^2Rsmaii-^i^{n,z,). 
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Subtracting these two equations from each other shows that the range A0 for drain-hole 
particles to reach Rsmaii is 

Lsin6' 
The solid angle of incident particles lying between 6 and 9 + /S.9 is 27r sin 9/\9 and so the 
fraction f^h of all incident particles with angular momentum L that become drain-hole 
particles and fall to ^ < Rsmaii is 



Uk - 4^ - -^ • (90) 



6.6. Drain hole dynamo power 

The strength of the equilibrium radial electric field produced by drain-hole particles 
can be estimated as follows: Before any drain-hole particles accumulate at small f , there 
is no radial electric field, but as the drain-hole particles accumulate, the radial outward 
electric field will develop. The force due to the radial outward electric field will oppose 
the gravitational and magnetic forces causing the inward motion of the drain-hole particles. 
The balance between these opposing forces is quantified by the radial equation of motion. 
In cylindrical un-normalized coordinates the radial equation of motion governing drain-hole 
particles with f ^ = in the z = plane (where B^ = due to 2;-antisymmetry of /(r, z)) is 

rud (f - r02) = g, (^Er + r^B,^ - ^^^^^. (91) 

i?2 is approximately uniform at small r so ^ ~ irr'^Bz at small r in which case the drain-hole 
particle condition P^ = mdr'^(f) + qaip/^TT = implies <p = —QdBz/'^^^d- On eliminating in 
EqjnH the radial equation of motion governing drain-hole particles is 

r - Id F r-Jd£l MG 

rrii Aval ^ f92l 

~ S<LE -MG ^ ' 

where the second line is for small r. When Er = 0, the drain-hole particles fall inwards with 
gravitational acceleration, but as Er builds up because of accumulation at small r of fallen-in 
positively charged drain-hole particles, Eq|92] shows that this electric field will oppose the 
gravitational force and retard the infall. The drain-hole particles will continue to fall in and 
accumulate, thereby increasing Er until the radially outward repulsive electrostatic force due 
to Er becomes so strong as to balance gravity and cause f to vanish. Thus, gravitational 
force is balanced by the radially outward force from the space charge electric field of the 
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accumulated positively charged drain-hole particles. The saturation electric field is 

Using Er = —dV/dr, integration of Eql93] from large r to the jet radius Vjet gives the 
voltage at the jet to be 

v.. = '^^i^. (94) 

Qdrjet 

The jet electric current corresponds to the charge per second carried inward by the 
drain-hole particles. The number of drain-hole particles accreting per second is Mdh/fna 
where M^h is the mass accretion rate per second of drain-hole particles and rrid is the mass 
of an individual drain-hole particle. Thus, the poloidal electric current is 



Ijet = qdMdh/rrid. (95) 



The jet electric power is 



MdhMG , , 

P,et = hetV.et = —^ (96) 

rjet 

which is just the rate at which gravitational potential energy is released by drain-hole parti- 
cles falling from large r to the jet radius. The drain-hole dynamo converts the gravitational 
energy released from accretion into electrical power suitable for driving bipolar jets that 
are moving away from the z = plane. The jet power is proportional to both the central 
object mass M and the drain-hole mass accretion rate Mdh- Paper I showed that the dust 
mass accretion rate can be a substantial fraction of the total mass accretion so Pjet can be a 
substantial fraction of the power of all accreting material. The jet power accelerates the jet 
material to escape velocity and so is equal to the power available from accreting drain-hole 
dust grains. Thus, assuming that the axial starting point for jet particles is of the order of 
rjet, the power required to drive the jet particles to escape velocity is Pjet = MjetMG/rjet 
and so the jet mass outfiow would be approximately equivalent to the drain hole particle 
accretion rate, i.e., Mjet — Mdh- The particles in the jet would not, in general, be the drain 
hole particles, but instead would be plasma magnetohydrodynamically accelerated using the 
drain-hole accretion as the power source. Assuming L ~ 1 and Rsmaii = fjet/a ~ 0.1 in 
Eqj90| the fraction of retrograde particles that are drain hole and able to reach Vjet would be 
fdh = 0.4; the fraction of combined retrograde and prograde dust grains would thus be 0.2. 
The example in paper I showed that because of differences in proportional slowing down, 
the dust accretion rate would be enriched to be 20% of the gas accretion rate. This gives 
Mdh/Mg ~ 0.2 X 0.2 = 0.04 and so predicts a jet power that would be about 1/25 of the power 
associated with all accreting dust and gas. This ratio of outfiow power to accretion power is 
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consistent with the estimate given by iBacciotti et al. I (120041 ) using HST observations of T 
Tauri jets. 



7. Torque and angular momentum 

An important question repeatedly asked about accretion disks and jets is the role played 
by jets in satisfying conservation of mechanical angular momentum of the accreting material. 
It will now be shown that mechanical angular momentum is exactly conserved in our model. 

Because of axisymmetry, the canonical angular momentum of the j^^ charged dust grain 

Ptd,j = ^+^^1 + Ze^i^, z)/27r (97) 
and the canonical angular momentum of the k^^ electron 

PI, = m,rvl-e^P{r,z) /27r (98) 

at any position r,z are both invariants, i.e., Pt^ = const, and Pf = const. In general, the 
dust grains and the electrons at any position r,z will have quite different values of v^ but, 
in order for the plasma to be macroscopically quasi-neutral, there must be approximately 
Z electrons adjacent to each dust grain. Because electron and dust grain trajectories differ, 
these neighboring electrons will typically not be the original sibling electrons photo-emitted 
when the dust grain became charged. 

The radial and axial velocities of a specific dust grain or electron can be written as 
f^'-' = dr^'^ /dt and f^'-' = dz"'^ /dt where r"'^{t) and z"'^{t) are the position of the j^^ 
particle of species a. Since P^ is conserved for each individual particle, the time derivatives 
of the P(^'s of a dust grain at a location r, z and its neighboring neutralizing Z electrons 
respectively give dP^^Jdt = and dP^i^/dt = 0. Using dip/dt = Vrdip/dr + Vzd4'/dz for the 
time derivative of ip measured in the particle frame, respective time derivatives of EqsJ971 
and [98] give 



'' 'Zm^rvf) = 1^ %'/ + %f ] . (100) 



and 

- [Zm^rv^ j - — y—v, + — t. 

Using Bj. = ~{27Tr)^^dtp/dz and B^ = {2'Kr)^^dil)/dr from Eql7]and summing Eqs 
and IIOOI over the dust grains and their associated Z neutralizing electrons at location r, z 
gives 

—^ = r {J,Br - JrB,) = r0 ■ JpoixBpoi (101) 
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where Jr-, Jz are the respective radial and axial current densities and L^ is the total me- 
chanical angular momentum density taking into account both dust grains and electrons. 
Thus, from the macroscopic point of view there is a torque about the z axis, namely 
z-r X F = z X [rf + zz) ■ {3poi x 'Qpoi) = r(f)-3poix'Bpoi acting to change the local mechanical 
angular momentum density. 

On the other hand, using Eqs|7] and [5] it is seen that when this torque is integrated over 
the entire volume to infinity, 

/ dh r0 ■ Jpoi X Bpoi = J d^r r^V(j) ■ (^ VJ x V0) x (^ V^ x V0) 

= 4^ / d^rr'^Vcj) x ( VJ x V0) ■ ( VV' x V0) 

= J^/ rf3rV/-(V^x V0) (102) 

= ^/ rf3^V-(/(VV^xV0)) 

= 

since both / and V^ vanish at oo. Thus, the total mechanical angular momentum of the 
system is exactly conserved because there is no net torque applied to the whole system. 

/ is an odd function of z and ip is an even function of z, and in the jet Jpo; is nearly 
parallel to Bpoi- This suggests the following generic form for the current 

/xo/(r,z) ~ AV'(r,z)tanh ( — — j (103) 

where h{r) represents the height of the accretion disk at radius r and A, the current per 
flux, has units of inverse length. Thus far from the z = plane, Eq J103l has the form 
/io^/V' = Asign(2;) so that the jet above the z = plane has the opposite handedness of the 
jet below the z = plane. The parameter A is closely related to the current per flux in a 
force-free system (i.e., a system satisfying V x B =AB), but differs slightly because here A 
refers to just the ratio of the poloidal current to the poloidal flux. Using Eq JlOll it is seen 
that the density of MHD torque about the z axis is of the generic form 



dt Att'^t V dr dz dz dr / 



X 



ifiOTT^r 



{-I: [^(r, .) tanh (^y)] |^ + fe [^(r, .) tanh (^y)] f } 



(104) 



We assume that dh/dr -C 1 so the radial scale length at which h changes is much larger 
than h. Also, from symmetry dip/dz = on the 2; = midplane. Together, these conditions 
imply that near the midplane the last term in Eq J104l dominates so near the midplane 

dL^ _ A V'(r, z) dj) ^ Xi){r;z) Bz{r,z) 

dt 4/iovr2/ir cosh^(2;//i) (9r 2TCjj,Qh cosli^iz/h)' 
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As seen from Eq J102l the torque density is proportional to — V/ x V(j) ■ Vip ~ — "rJpoi ■ V^ 
and so is positive for poloidal current flow away from the poloidal field magnetic axis and 
negative for poloidal current flow towards the poloidal field magnetic axis. The torque density 
vanishes as r — >■ and as r — > cxd and at the poloidal field magnetic axis because r VV' — *■ 
at r = 0, Jpoi — i> at r = oo, and Vz/' = at the poloidal field magnetic axis. The direction 
of poloidal current flow is shown in FigJTSl This torque acts on the drain hole particles 
and their associated electrons since these particles are the carriers of the poloidal current as 
sketched in Fig Upland FigIT51 Unlike the drain-hole particles, no torque r0- J x B about the 
z axis acts on the Speiser particles because the current associated with the Speiser particles 
is in the (h direction. 



8. Conclusions 

We have shown that charging of coUisionless dust grains incident upon a star causes 
the dust grain orbital dynamics to change from a relatively simple Kepler form to more 
complicated motion involving competition between magnetic and gravitational forces. This 
competition gives rise to five qualitatively different types of orbits. Two of these, the retro- 
grade and prograde cometary orbits are just perturbations of Kepler cometary orbits. The 
orbit of a particle where magnetic forces overwhelm gravitational forces is just a Larmor (cy- 
clotron) orbit and in this case the particle is constrained to remain within a poloidal Larmor 
radius of a poloidal flux surface in a manner similar to tokamak confinement. Particles where 
magnetic and gravitational forces are comparable can have two very different types of orbit 
depending on whether the incident particle is prograde or retrograde. Prograde particles 
of this latter type develop Speiser orbits; these orbits are paramagnetic with respect to the 
poloidal magnetic field and so can be the source of the poloidal magnetic field. Retrograde 
particles having comparable magnetic and gravitational forces can have a peculiar behavior 
whereby centrifugal force is eliminated with the result that the charged particle falls in to- 
wards the star along a spiral orbit. The accumulation of these "drain-hole" particles near the 
star provides a radial electric field oriented so as to drive the poloidal currents and toroidal 
magnetic fields of an astrophysical jet. 

This paper showed the existence of these different types of orbits, how their orienta- 
tion is suitable for generating the poloidal and toroidal magnetic fields associated with an 
accretion disk and astrophysical jet, and how questions of angular momentum conservation 
are inherently resolved. A future paper will investigate the quantitative values of dust grain 
parameters required to produce toroidal and poloidal fields in the accretion disk of a young 
stellar object. 
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Finally, we offer some remarks regarding the effect of deviations from axisymmetry. 
Tlie model presented liere assumed perfect magnetic field axisymmetry field whereas actual 
accretion disks are observed to have varying amounts of non-axisymmetry. This situation 
is analogous to toroidal magnetic fusion devices such as tokamaks, reversed field pinches, 
and spheromaks all of which are modeled to first approximation as being axisymmetric, 
but in reality have deviations from axisymmetry due to waves, turbulence, instability, and 
errors in machine construction. It is known from these devices that a modest breaking of 
symmetry does not invalidate the results of the axisymmetric model, but rather weakens 
the conclusions, e.g., instead of cyclotron-orbiting particles being perfectly confined to the 
vicinity of a poloidal flux surface, when there is deviation from axisymmetry cyclotron- 
orbiting particles can slowly wander away from the poloidal flux surface they started on. 
One would expect that deviations from axisymmetry in accretion disks would cause a similar 
transport of cyclotron particles across poloidal flux surfaces. Because symmetry breaking 
causes the canonical angular momentum of particles to change, it could be considered as 
being somewhat like a collision that changes the canonical angular momentum of each of two 
particles involved in a collision while conserving the total canonical angular momenta. Hence, 
deviations from axisymmetry would cause a jiggling of the canonical angular momenta of 
individual particles so that particles on the borderline between being drain-hole and cyclotron 
or on the borderline between being Speiser and cyclotron might spend part of the time (i.e., 
between jiggles) being one type and part of the time being the neighboring type. Similarly 
cyclotron particles that are on the borderline between being mirror-trapped and not mirror- 
trapped would, as they get kicked into and out of the mirror loss-cone, spend part of the time 
being mirror-trapped and part of the time not being mirror-trapped. However, at any given 
time there would be a certain fraction of particles of each type, i.e., a certain fraction would 
be cyclotron, a certain fraction would be drain-hole, a certain fraction would be Speiser, and 
a certain fraction would be cometary. 
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A. Derivation of generic poloidal flux function 

If all the toroidal current T^ is concentrated at the poloidal location r = Rq and z = 0, 
then the toroidal current density is 3 tor = 4>^(i,S{z)6{r — Rq). On defining 

^' - '''::. . (Ai) 



(Ro + ry + z^ 



analytic solution of EqE] using J tor = 4>^(i>S{z)6{r — Rq) gives (1Jacksonlll999l ) 



ij{r,z)= ^^/R:^[{2-e)K{k)-2E{k)] (A2) 

where E and K are complete elliptic integrals. Equation IA2I describes the situation where 
all the current density is concentrated at r = Rq, z = 0, i.e., the current fiows in a wire of 
zero cross-section located at r = Rq, z = 0. This equation can also be used to (i) describe the 
field observed at locations far from the poloidal field magnetic axis of a distributed current 
localized in the vicinity of the poloidal field magnetic axis and (ii) as the Green's function for 
a distributed toroidal current. This is because for an observer who is far from r = Rq, z = 0, 
the field of a distributed toroidal current localized near r = Rq, z = is indistinguishable 
from the field of a zero cross-section wire carrying the same total current. Equation IA2I has 
a logarithmic singularity at the wire location because the wire has infinitesimal diameter. 

Two analytic limits are of interest for Eq JA2[ The first is where r <^ Rq so 

, 2 4:RQr 



2 , ^2 



and the second is where r ^ Rq so 



m + z- 



2 4i?or 



(A3) 



k' ^ ^rt-2- (A4) 

The former gives the field near the loop axis and the latter gives the field at locations far 
from the current loop. In both cases k"^ is small compared to unity and so the small argument 
asymptotic expansions of the complete elliptic integrals can be used, namely. 



so for r <^ Rn 



and for r ^ Rq 
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Thus for small k, it is seen that (2 — /c^) K{k) — 2E{k) ~ Tck^/IG in which case 

^(^,^)-^7 ^m^. (A6) 

lim^(r,z)^ ""^°^^ ^°^\,^ (A7) 

lim^(r,z) ^ ^/^oX^ ^°''\;. - (A8) 

r>a^^ ' ^ 2 (^2 + ^2)3/2 V 7 

For purposes of discussion and also numerical computation, it is convenient to choose Eq lA6l 
to represent the poloidal flux of a generic toroidal current everywhere. Making this choice 
for the poloidal flux function (instead of the prescription given by Eq JA2p means that ilj{r, z) 
has a smooth hill-top at r = 2Rq rather than a logarithmic singularity at r = i?o and has 
the same behavior far from r = Rq, z = Q as does Eq ]A2[ 

Thus, a useful analytic representation for a nonsingular, physically realizable flux func- 
tion is obtained by recasting Eq JA6l in the form 

,/ X 27 (r/a) , ... 

V(r,2;)= — ^^^ — -V^o. (A9) 



8((^ + i)'+(i)^ 



3/2 



This has a maximum of ipQ aX r = a, scales as r^ for small r, and scales as r~^ for large r. 
Equation [6] can be used to calculate the associated toroidal current density which will be 
sharply peaked near r = a and 2; = 0. The ip{r., z) prescribed by Eq JAQl has the features 
that it provides a dipole-like field far from the z axis and a nearly uniform axial field near 
the z axis, corresponds to a realistic distributed toroidal current, has a well-defined poloidal 
field magnetic axis, is analytically tractable, and is convenient for numerical computation of 
representative particle orbits in a physically relevant magneto-gravitational field. 



B. Current associated with flux function 

Using Ampere's law to relate the toroidal current and the poloidal magnetic field it is 
seen that 

X^ = — i Bp,i-d\ (Bl) 

where the contour C links the total toroidal current X^. By letting the line integral go to 
infinity in the radial and z directions it is seen that only the portion of the line integral along 
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the z axis makes a finite contribution so 






li-^-^o^iU^X 



167rAfo ^^^^^'"-*0 rar I a2 J-oo "777 ix2 /^n2\3/2 



Defining h = r/a + 1/2 and z/a = fesinh-i? the z integral can be expressed as 

TOO dz _ „ r°° fecoshi?di9 

J-OO //. . -,\2 /,x2\3/2 — " J_ 



= f[tanh^]!^^ (B3) 



2a 

(r/a+1/2)^ • 



Since 



lim -— - 

r— >o r (7r 



r^ 2a 

a^ (r/a +1/2)^ 



16 
a 



(B4) 



the total toroidal current is 

I, = ^^. (B5) 

vra/io 



C. Review: Distinction between diamagnetic (adiabatic) orbits and 

paramagnetic (Speiser) orbits 

C.l. Diamagnetism of cyclotron (Larmor) orbits 

We first review charged particle motion in a uniform magnetic field B = B^z (so ip = 
BzTir'^) and no electric field; orbital motion in more complex fields will be discussed later. 
The particle motion is prescribed by the Lorentz equation 

m^— = q^v X B^z . (CI) 

If the particle is restricted to the z = plane, the respective radial and azimuthal components 
of EqEUare 

TTig. ( f — r0^ ) = q^-rcpBz (C2) 



We consider circular motion (i.e., cyclotron or Larmor orbits) so r = const, in which case 
Eq JC3l gives (p = const, and Eq JC2l then becomes 

= -OOca (C4) 
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where ujca = qaBz/rricr is the signed cyclotron frequency. The minus sign in Eq lC4l indicates 
that cyclotron motion is diamagnetic. Thus if a gyrating charged particle is considered 
as a 0-directed current, the polarity of this current is such as to create a magnetic field 
which opposes the initial field B^, i.e., cyclotron orbits tend to depress the value of ijj. The 
diamagnetism of cyclotron orbits means that cyclotron orbits cannot be the source for the 
assumed poloidal magnetic field ipi"^, z) nor the means by which this field is sustained against 
dissipation. 



C.2. Adiabatic orbits 



When the magnetic field is non-uniform or there are electric fields, and if these addi- 
tional features are sufficiently weak that to lowest order the Larmor orbit (cyclotron orbit) 
description is approximately correct, then additional charged particle motions occur which 
are superimposed on the Larmor orbits v/^(t); these additional motions are adiabatic in the 
sense of classical mechanics. Defining f y as the velocity component parallel to the magnetic 
field and v^ as the component perpendicu lar to the magnetic field these motions are the 
standard drifts (JLongmird 119671 : IChenlll984j ). namely the E x B drift v^; = E x B/i?^, the 
polarization drift Vp = m„q~^B~'^dFjj_/dt, the curvature drift Vc = —m^v'uB ■ VB x B/q^B"^, 
and the grad B drift vys = —fiVB x B/q^B"^ where /i = mrjv']_/2B is the magnetic moment, 
an adiabatic invariant. There is also a 'force' drift v^t^ = F x B/q^B"^ where F is a generic 

1/9 

non-electromagnetic force, which here is gravity, so F = mMCV {r'^ + z'^) . Taking into 
account all these drifts, the velocity of an adiabatic-orbit charged particle becomes 



V\lB + YL^{t) +^ + 



dE 






+ 



B'^ ' q^B^ dt 
m„MG'^ ( 1 



rB^ 



\Jr'^-\-z 



mCTfjjB-VBxB 

X B. 



(C5) 



The last four drifts in Eq ICSI explicitly involve q^ and thus produce macroscopic currents. 
When these currents are summed and, in addition, diamagnetic current is taken into account, 
the result is equivalent to th e MHD equation of mot i on where the pola rization drift plays 
the role of the inertial term (JGoldston fc Rutherford! Il995l : iBellanI |2006| ) . The ideal MHD 
concept of frozen-in flux is directly equivalent to \x conservation because /i conservation 
corresponds to conservation of the magnetic flux linked by a cyclotron orbit. Thus, the ideal 
MHD concept of frozen-in flux is based on the adiabatic invariance of cyclotron orbits. 

The poloidal flux function specified by Eqj9] corresponds to a magnetic field generated 
by a toroidal current flowing in the positive direction (counterclockwise direction); the B^ 
component of this field is positive for r < a and negative for r > a where a is the location of 
the poloidal field magnetic axis. The poloidal magnetic field has both curvature and gradients 
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so that away from field nulls, particles should have parallel motion and cyclotron orbits 
together with superimposed curvature and grad B drifts. Figure [H] shows the numerically 
calculated orbit of a particle located in the z = plane in a magnetic field prescribed by 
Eqini and located inside the poloidal field magnetic axis (indicated by dashed circle). It is 
seen that the particle makes cyclotron orbits with a superimposed drift due to curvature and 
V-B. The cyclotron orbit is clockwise consistent with the assertion that cyclotron motion 
is diamagnetic. Figure [15] shows the situation for a particle located at a radius outside the 
poloidal field magnetic axis. The sense of the cyclotron orbit is now reversed as is the polarity 
of Bz so the cyclotron orbit is again diamagnetic. For both inside and outside particles the 
drift motion is clockwise and so opposes the original toroidal current creating the poloidal 
flux and so the curvature and V-B drift motion can also be considered diamagnetic. 

The current associated with the gravitational force drift is 

J. = E «<,,.v, = ^V (;p^) X B (C6) 

where g = Yl ''^0^^ is the mass density. 

From a macroscopic (i.e., MHD) point of view, the force associated with the gravita- 
tional drift current exactly balances the gravitational force component perpendicular to the 
magnetic field since 



J.xB = ^(v(^;^)xbUb 



-qMGVa 



1 



\/t'^+z 



(C7) 



If / = on the z = plane (as is consistent with astrophysical jet symmetry used by iLovelace 



( 1l976l )). the gravitational drift is not defined on the poloidal field magnetic axis because B 
vanishes on the poloidal field magnetic axis and the theory of particle drifts fails. In other 
words, going from the first to the second line in Eq JC7l at the poloidal field magnetic axis 
would involve dividing zero by zero (since 5 = on the poloidal field magnetic axis). 

When summed over species, the curvature and grad B drifts correspond to currents 
which balance macroscopic pressure gradients (when diamagnetic current is included) and 
the polarization current correspon ds to the ine rtial term in the MHD equation of motion. 



This analysis shows, as discussed in iBellanI ( 120071 ) . that plasma particles undergoing cyclotron 



motion and drifts do not have Keplerian orbits. It also shows that the poloidal field magnetic 
axis is a special place where conventional particle drift theory fails. 
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Fig. 14. — Orbit of a positive particle in the z = plane located inside the poloidal field 
magnetic axis (indicated by dashed circle), coordinates are normalized to the poloidal field 
magnetic axis radius. Bz is positive inside the circle and negative outside. 
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Fig. 15. — Orbit for a positively charged particle located in the z plane outside the poloidal 
field magnetic axis (dashed circle). 
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C.3. Non-adiabatic motion: the Speiser orbit 

An extreme form of magnetic non-uniformity occurs wliere tlie magnetic field reverses 
direction. In tliis case an orbit quite distinct from the cyclotron orbit and its associa ted adi- 



abatic drifts occurs. This non-adiabatic orbit, called a meandering or Speiser orbit (ISpeiser 



19651 ). consists of semi-circles of counterclockwise motion interspersed with semi-circles of 



clockwise motion. 

A numerically calculated Speiser orbit for a positively charged particle in the z = plane 
is shown in FigJT6l The particle oscillates across the poloidal field magnetic axis between the 
inside region where Bz > and the outside region where Bz < 0. The result is a net coun- 
terclockwise motion so, in contrast to cyclotron orbits, Speiser orbits are paramagnetic. The 
paramagnetism of Speiser orbits has been considered an important aspect of current sheets 



in Earth's magnetotail, [e.g., see IZelenyi et al. I ( |2000| )]. but to the author's knowledge this 



paramagnetism has not been previously considered in the axisymmetric three-dimensional 
geometry discussed here which is relevant to accretion disks and astrophysical jets. In par- 
ticular, we will show that the poloidal flux function can be considered as a consequence 
of Speiser orbits such as shown in FiglTBl Speiser orbits are not consistent with the drift 
approximation (i.e., E x B drift, grad B drift, curvature drift, etc.) because the drift 
approximation is based on the assumption that, to lowest order, the particle is undergoing 
cyclotron motion. The inconsistency between Speiser orbits and the drift approximation 
is obvious when one considers that the drift approximation fails where B reverses polarity 
whereas Speiser orbits depend on this reversal. 

If motion in the z direction is also allowed, then because B^ also reverses at the poloidal 
field magnetic axis, the particle can also oscillate vertically across the poloidal field magnetic 
axis to make vertical Speiser orbits. The combined r and z Speiser motion means that 
particles moving at an arbitrary angle across the poloidal field magnetic axis will reflect 
from interior surfaces of the nested poloidal flux surfaces concentric with the poloidal field 
magnetic axis. These nested poloidal flux surfaces can thus be imagined as the walls of 
a toroidal tunnel and the Speiser orbit particles can be considered as reflecting from the 
interior walls of this toroidal tunnel while moving in the counterclockwise direction to trace 
out paramagnetic orbits and create poloidal flux. 



D. Equation of motion and its solutions 

The Hamiltonian orbit analysis presented here shows that photo-emission creates new 
effective potential barriers. The topography of these barriers depends on a combination of 
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Fig. 16. — Speiser orbit. The charged particle bounces back and forth across the field null 
at the poloidal field magnetic axis resulting in a counterclockwise (i.e., paramagnetic) orbit. 
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environmental factors, particle properties, and the location of the charging. Representa- 
tive orbits obtained by numerically integrating the equation of motion have been presented 
and are consistent with the predictions of the Hamiltonian theory. We outline here the 
derivation of the dimensionless equation of motion; this derivation gives insights into several 
fundamental issues regarding the dynamics, especially the influence of initial conditions. 

The equation of motion for a charged particle in a combined electromagnetic and grav- 
itational field is 

ma-^ = qA^ + ^x^)+m<rMGV^^==. (Dl) 

Using Eq{T] for the magnetic field, the equation of motion can thus be written as 



+ l^f x(ixi + /^o/i) (D2) 



Then, using the definitions given in EqJSTl the equation of motion can be expressed in 
dimensionless form as 



where 






is the dimensionless electric field and 

V = ^^-^ (D5) 

is the dimensionless electrostatic potential. 

Equation ID3I clearly shows that the dynamics change from being gravitationally domi- 
nated to being magnetically dominated according to the ratio {oJca) /^o- The possibility of 
complex interactions between gravitational and magnetic forces when (cUco-) /^o is of order 
unity is also evident. The coefficient fiol/aii {B^) is related to the pitch of a twisted field. 
The Hamilton-Lagrange formalism shows that / plays a subservient role for particle orbits 
compared to ip because canonical angular momentum depends on ip, not /. However, large I 
increases \B\ and so contributes to the effective potential /x|-B| thereby providing additional 
possibilities for localization. Thus, if /i is large, a particle is not only constrained to stay on 
a constant ip surface, but is additionally constrained to stay out of regions on this surface 
where fi\B\ is large. If poloidal currents fiow, then the associated Jpo^xBtor forces drive jets 
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which inflate and distend the ip surfaces. Tlius, tlie orbits will depend indirectly on / when 
the jet dynamics alter the shape of the constant ip surfaces. 



Using the relations 



d4i </> —xy+yx 

dj) _ dii_^ _|_ d^^ _ d^ ( xx+yy \ , dij_^ ^ > 

(9x df dz df \ /^iipp' j dz ' 



the normalized equation of motion can be expressed in Cartesian coordinates as 



d?5L _ dV 



( dtp i, dtp r 
df dzf I . - 

. ^qJ i-xy+yx) I {^') 

jr+zz) 

|^2+52|3/2 



where r = xx + yy and r = ^/x"^ + y"^. Equation ID7I is in a form suitable for numerical 
computation and has been used to provide the orbital plots shown earlier. 

At this point it is convenient to use the generic poloidal flux function given by Eqj9] so 
the unity-maximum, dipole-like, normalized flux function will be 



27f' 



with 



v^v ?^; 


8((f+if + .-f'' 


dip 


27f(f + 1 + 4^2- 2f2) 


df 


i6((.-^ir + rf'= 


dip 


Slf^z 



and 

Thus, the normalized poloidal magnetic fleld components are 



(D8) 



(D9) 



(DIO) 



81rz ,^ 

Br = -TT^ (Dll) 
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We now consider the problem of establishing appropriate initial conditions for an in- 
coming neutral particle. For purposes of starting a computation we assume the particle is 
located at some initial radial position pQ in the orbital plane such that po > Ppericenter where 
Ppericenter IS givcu by Eq|271 Solviug Eql^S] for the initial inward radial velocity gives 



_ 2i?-^ + -. (D13) 

*' Po Po 

and the corresponding initial orbital frame azimuthal velocity is 

f;,o = -. (D14) 

Po 

Equation [30] can be solved for the initial polar angle in the orbital frame as 

V = a + cos-i ( ]-^'Ip^ \ . (D15) 

Using Eq|3T] the initial orbital frame Cartesian coordinates are thus 

J.I _ L'^ cos r)o 



1-^1+2L2 Hcos(rio-a) 

y' = ^_''^''^_^» (D16) 



~/ 



Z' 



= 0. 



1 - \/i+2L2 H cos (-no-a) 



The orbital frame Cartesian velocity components are related to the orbital frame cylindrical 
velocity components by 

Vx'o = Vpo COS ?7o - Vr^Q siu ?7o 

Vy'o = Vpo sin r]o + v^o cos r]o (D17) 

Vz'O = 0. 

We now take into account that the orbital frame Cartesian coordinate system is rotated by 
the angle of inclination 6 about the x axis with respect to the lab frame coordinate system. 
The X and x' components of both position and velocity are the same in the two frames but 
the y and z components are related by 

y = -z' sine + y' cose 
z = z' cose + y' sine. 

Since z' is by definition zero in the orbital frame, the initial lab frame Cartesian coordinates 
are 

Xq = Xq 

yo= y'ocose (di9) 

zq = 2/0 sin e. 
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Since v'.^ is similarly zero in the orbital frame, in analogy to Eq JDlQl the initial lab frame 
Cartesian velocities are 

VxQ = Vx'O 

VyO = VylQCO?,6 (D20) 

VzQ= VyiQsinO. 

Thus, if one wishes to start the numerical computation at the radius po on the trajectory 
of an incoming particle with orbit parameters {H,L, 6, a}, Eqs JD13l ID14t [PTSt ID19I and lD20l 
give the appropriate initial position and velocity lab frame Cartesian components. Before 
charging, the orbits are degenerate with respect to choice of ^ or a, but after charging there 
is a strong dependence on these two angles. In particular, if < 6' < 90° the orbit is 
prograde and Speiser type orbits are possible if the charging occurs near the poloidal field 
magnetic axis. On the other hand if 90° < < 180° the orbit is retrograde and drain- 
hole orbits are possible. Thus, a subclass of prograde incident neutral particles transform 
upon charging into the toroidal-current /poloidal-field dynamo while a subclass of retrograde 
neutral particles transform upon charging into the poloidal-current/toroidal-field dynamo 
that drives a bipolar astrophysical jet. Because 9 and a also affect the angle between the 
velocity vector and the magnetic field at charging, 9 and a affect the value of p and hence the 
extent to which accreted particles with cyclotron orbits will be mirror trapped to subregions 
of constant ip surfaces. For example, if a = then variation of the angle of inclination 9 
for a given Ppericenter, and charging radius R^: will determine whether the charged particles 
created upon disintegration of an incoming neutral particle will be normal trapped particles, 
untrapped particles, drain-hole particles, or Speiser particles. 
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